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Paul Scherrer 7ojahrig 


Am 3. Februar 1960 begeht Prof. Dr. PAuL SCHERRER, Direktor des Phy- 
sikalischen Institutes der ETH, seinen 70. Geburtstag. Dieses Ereignis bedeutet 
gemass den tiblichen Vorschriften einer Hochschule mehr als nur einen Markstein 
im physischen Leben eines Menschen, sondern bedingt meistens den Riicktritt 
von Amt und Wiirden und im Falle SCHERRER das Ende einer mehr als 30jah- 
rigen Epoche der physikalischen Lehr- und Forschungstatigkeit, welcher der 
Jubilar sein ureigenstes Geprage aufdriickte. 

1920 kam SCHERRER zusammen mit dem grossen Gelehrten DEBYE von 
Géttingen nach Ziirich, um hier die zweite Lehrstelle fiir Experimentalphysik 
an der ETH zu bekleiden. Es folgte eine selten gliickliche Zeit am Physikalischen 
Institut der ETH, die getragen war von den beiden Persénlichkeiten DEBYE 
und SCHERRER. Das lebendige Zusammenspiel der beiden, in ihren Grundziigen 
so verschiedenen, Gelehrten liess seine befruchtende Wirkung weit in die Welt 
ausstrahlen und Ziirich in kiirzester Zeit zu einem Physikzentrum internatio- 
nalen Formats heranwachsen. 

Leider fand diese herrliche Episode schon 1927 ein Ende, indem DEBYE 
sich entschloss, dem an ihn ergangenen Ruf an die Universitat Leipzig Folge 
zu leisten. An der ETH wurden daraufhin die beiden Hauptprofessuren fiir 

_ Experimentalphysik in eine zusammengelegt und diese SCHERRER tibertragen, 
dessen Lehrerfolg sich bereits damals eines grossen Rufes erfreute. Damit war. 

_ der Start fiir die weitere Entwicklung von Forschung und Unterricht auf dem 
Gebiete der Physik vorgezeichnet, und zwar in einer Eigenart, wie sie nur 
einmalig anzutreffen ist. Wenn der Schreibende mit besonderer Betonung auf 

_ die Jahre hinwies, wo das Doppelgestirn DEBYE—SCHERRER das physikalische 
Geschehen in Ziirich beherrschte, so auch deshalb, weil es ihm personlich ver- 
génnt war, diese Zeitspanne als Assistent von DEBYE mitzuerleben und sie ihm 
wie ein leuchtendes Kleinod unausléschlich in Erinnerung blieb. 

Die nun folgenden drei Dezennien zeugen von einer unerwarteten, machtigen 

Entwicklung der physikalischen Disziplinen an der ETH. Sicherlich hat auch 
WOLFGANG Pautt, der sich dem Jubilar kurze Zeit nach seiner Ernennung 
zum Leiter des Physikalischen Institutes der ETH, als Ordinarius ftir theore- 

_ tische Physik, beigesellte, entscheidenden Einfluss auf das physikalische Leben 

in Ziirich; im engeren Sektor der Experimentalphysik war es aber SCHERRER 

-allein, der dem Geschehen seine ihm charakteristischen Ziige aufpragte, alles 

iiberstrahlt und beschattet von einer beispiellosen Meisterschaft in der 
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Darstellungskunst, wie sie sich vor allem in den Vorlesungen SCHERRERS in 
seltener Bliite offenbaren konnte. Jede neue Erkenntnis auf dem Gebiete der 
Physik wurde unmittelbar nach ihrem Bekanntwerden in das Pensum der Vor- 
lesungen iiber Experimentalphysik eingebaut und in einer begliickenden, leben- 
digen und einfachen Weise vorgetragen, die die Zuhérerschaft stets in helle 
Begeisterung versetzte. Eine Begeisterung, die mitgeholfen hat, die Kunde 
von der Meisterschaft des Jubilars weit tiber die ganze Welt zu verbreiten. 

Die Wiirdigung SCHERRERS als Forscher und Gelehrter findet ihren schén- 
sten Ausdruck in einer im Birkhauser Verlag, Basel, erschienenen Festschrift 
zur Feier seines 70. Geburtstages, in welcher sich zahlreiche seiner Schiiler und ~ 
Mitarbeiter zum Worte melden und Zeugnis ablegen von der Vielgestaltigkeit 
seiner wissenschaftlichen Tatigkeit. Hier sei nur darauf hingewiesen, dass diese 
Tatigkeit namentlich Fragen des strukturellen Aufbaues der festen K6rper 
und der Deutung ihrer physikalischen Eigenschaften beriihrte und sich auch 
auf eine Fiille von Problemen aus dem Gebiete der Kernphysik erstreckte, 
wobei SCHERRER schon frithzeitig, unmittelbar nach Bekanntwerden der ersten 
Befunde und Entdeckungen, die fundamentale Bedeutung dieses neuen For- 
schungszweiges in seiner vollen Tragweite erkannte. 

In tiberzeugende Weise hatte SCHERRER schon kurz nach dem Zweiten Welt- — 
krieg begonnen, in Vortrégen vor Kreisen der schweizerischen Wirtschaft und — 
Industrie wie auch der breiteren Offentlichkeit auf die iiberragende Rolle auf- 
merksam zu machen, welche die Kernenergie als Energielieferant der Zukunft — 
spielt. Es ist weitgehend sein Verdienst, wenn heute in der Schweiz eine Kern- © 
reaktoranlage in Betrieb steht, die sowohl in bezug auf die Energiegewinnung 7 
und Isotopenverwertung als auch fiir die Heranbildung eines Stabes qualifi- — 
zierter Fachleute bereits zukunftsverbundene Arbeit leistet, zu Nutz und 
Frommen der schweizerischen Wirtschaft. i 

In seiner Eigenschaft als Redaktor der ZAMP verbleibt dem Schreibenden — 
noch die angenehme Pflicht, dem Jubilar fiir seine bei der Griindung der — 
Zeitschrift vor 10 Jahren bekundete Bereitwilligkeit zum Beitritt in die Redak- 
tionskommission und fiir seine in der Folge getatigte Mithilfe wahrend eines 
ganzen Dezenniums herzlich zu danken. R. SANGER © 
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On the Stress Distributions Due to Force Nuclei in an Elastic 
Solid Bounded Internally by a Spherical Hollow and 
in an Elastic Sphere 


By Derek CoLtins, Newcastle upon Tyne, England!) 


1. Introduction 


Although solutions for the problem of a force nucleus acting at a point in 
a semi-infinite isotropic elastic solid with stress free surface have been known 
for some time [1-3]?), the corresponding problem of a force nucleus acting 
at a point in an infinite elastic solid containing a spherical hollow appears to 
have received little attention. The simplicity of the solutions of the problem 
for the semi-infinite solid and of similar problems in the Stokes flow of viscous 
fluid [4] suggest that the stresses in the problem for the solid containing a 
spherical hollow ought to be expressible in terms of elementary functions only. 
In an attempt to find such a solution I obtained general expressions for the 
stresses due to elastic distributions in an infinite solid containifig a spherical 
hollow in terms of the stresses due to the same distributions in an infinite 
unbounded solid, these being given in a recent paper [5], henceforth referred 
to as SH. If the elastic distribution is taken to be a force nucleus, expressions 
are obtained for the stresses, which, while in closed form, contain elementary 
functions and integrals which it does not appear possible to evaluate as ele-. 
mentary functions. These integrals could however be evaluated numerically 
to obtain the stresses on the boundary of the hollow. However, in practice, it 
proved more convenient to obtain expressions for the stresses as infinite series 
of spherical harmonics and use these for numerical work. 
In section 2 of the present paper I give expressions for the stresses due to 
_a force nucleus acting at a point of an infinite unbounded solid and in section 3 
expressions for the additional stresses which must be added to these to give the 
stress distribution due to the nucleus in a solid containing a spherical hollow, 
these expressions being infinite series of spherical harmonics. Expressions 
in closed form are also given for the stresses on the boundary of the hollow, the 
expressions in this form for the complete stress distribution being omitted be- 
‘cause of their length. Numerical examples are given to show the variation of 
the stresses over the boundary of the hollow. In section 4 I consider the 


1) University of Durham, King’s College, Department of Mathematics. 
2) Numbers in brackets refer to References, page 16. 


* 
F 


4 DEREK COLLINS ZAMP 


problem of an elastic sphere in equilibrium under the action of two equal and 
opposite force nuclei at internal points of the sphere, expressions for the stresses 
as series of sperical harmonics and numerical examples being given. 


2. Force Nucleus in an Infinite Solid 


We first give the stresses due to a force nucleus acting at a point of an 
infinite unbounded solid. If (x, y, z) are the cartesian and (7, 0, y) the spherical 
polar co-ordinates of a general point referred to an origin O, the positive z-axis 
being the line 9 = 0, we suppose the force nucleus of magnitude F acts at the 
point O,, cartesian co-ordinates (0, 0, f), f > 0, and is directed along the z-axis 
in the positive sense. LovE [6] then gives the cartesian components (up, Vp, Wo) 
of the displacement at the point (x, y, z) as 


ty = AD y= AE, ma (+ 35), en 


where R? = x* + y? + (z— f)?, F= 16a 4 A (1—~y), being Potsson’s ratio 
and yu the rigidity. ping iro 

We can now find the stresses 77%, 7 0, 0 0) and @ Gy in terms of elementary 
functions from the stress-strain relations. The notation for the stresses is that 
of SH, the suffix 0 indicating stresses in the unbounded solid. When 7 < f, 
these stresse?’can be expanded in series of spherical harmonics as 


fe n=0 (2 + 3) ia (2 2) 
m(n—1)(n+3—4n) /7r\n-2 f 
me ees ime sab | 

Oy = — 2HA Srp (Pt 2n—1+ 20) (rym | 
Fina 248 Ss [ot tan ots 2a) (1) 


00,— 2" bbe fata tint t an) (op 


eel aee as een |e 
(9 -t 345 v\n-2 
ena CG)” | Be (e088 ad 


2A coté nae (n+ 5 — 47) v\n 
cea & | gata (7) 
ant t SEA) Stee 
(2 — 1) (+ 
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n=1 
am (n + 3. —4 +) ge" 
(2% — 1) (5 
where P,(cos@) is the Legendre function of the first kind of degree » and 


T;(cos@) the associated Legendre function of the first kind of degree » and 
order 1, defined by [7] 


J "] TE (eos), 


Ti(o) =—(1 — 0212 al) g = cos. 


When r > f, the expressions for the stresses are obtained from equations (2.2) 
to (2.5) by replacing m by — (n + 1) in the m-th term of each series. 


3. Force Nucleus Acting at a Point in an Elastic Solid Containing an 
Unstressed Spherical Hole 


We now consider a solid bounded internally by an unstressed spherical 
hollow r = a (a x <A), centre O. To obtain expressions in closed form for the 


stresses 77, ry, 66 6, and @q@,, which must be added to the stresses 77», 7 r Oo, 


A6 6, and ¢ gy respectively to give the stress distribution due to a force nucleus. 
acting at the point (0, 0, f) of this solid, its axis directed radially outwards from 
the centre of the hollow, we use equations (3.16), (3.18), (4.5) and (4.7) of SH. 
In the last two of these equations and in equation (4.6) the expressions G(a) 


and H(a) should each be multiplied by a factor a~*. The stresses 77 and 7 6, and 


the stress combination 0) = 77) + 6 8, + Gy for the nucleus in the unbounded 
solid, which are to be substituted in these equations, are taken in the forms 
containing elementary functions and not as given by equations (2.2) to (2.5). 
The stresses obtained in this way are in closed form, containing elementary 
functions together with integrals which it does not appear possible to evaluate 
as elementary functions. Because of the length of these expressions, only 


those for the stresses 00 = 0 Oo + 0 6, and $9 = G%+ PV on the boundary 
of the hollow are given. These are 


66a) = #4 [FPO + 0,0) - 400) +76] G1) 
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and 

bi BY (i = Q,(0 
¢$(a) = cs we ee cos 6 — f) 4+ Q;(6) + 4 Q3(0) + 2908) , 32) 
where 


Rea ah ye ae 
The functions Q,(0); 7 =1,..., 6, are given by 
Q,(0) = (1 — 2 7) a® cosO + (1 — 2 9) a? f (1 — 4 cos?9) 


+ af?cos0 (—4+ 2 + (7-8) cos?6) + (4-24 — (5 — 40) cos?6), 


Q2(0) = (1 — 9) {2 (3 — 89) I (6-2, 5,1) 
= {1+ 12 4 — 16 7) s* Jle= 2,58} 
+ (3 — n) a {2 (25 — 64 7 + 40 7) Ic, s, 1) 
— (27 + 20 n — 128 7? + 80 7) s+! fic, s, 1)} 
+ fcosd {2 (— 51 + 183 » — 180 7? + 48 7) I(c — 1, s, 1) 
+ 3 (29-77 — 110 72 + 120 43 — 32 44) s4 J(e —1, 5, 1} 
+ (1 — 7) a® f’ cos6 {2 (1 — 4 7) I(c —1,s, 3) 
, — (5+2n—8 7%} s4 Jic—1,s, 3)} 
— (3 — n) a® f' cos6 {2 (5 —4 7) I(¢ +1, s, 3) 
+ (3-10 +87) s+ J(c + 1,8, 3)} 
+ 2 a8 I(c, s, 3) {47 — 70 4 + 16 7? — (27 — 44 1 + 8 97) cos?O} 
+ a* st J(c, s, 3) {3 — 102 y + 140 9? — 32 »8 
+ (9 + 56 4 — 88 4? + 16 73) cos? 6} | 


(0) = 2O— Ma 2289 ty (1 4. cos) — YF + Ry) (1+ cos). 


(1 — n) a? cos? 6 
2 | 


(7 + Ay) 


0 0 
i oe [cosec® : z ok 
2 R cose? 


ae ae [(R, — a) (a (ie n) Si ”) f cos6) 


— 7 (3 — n) a? cos? — (11 — 3») af’ cos0 + (17 —5y)f 
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Q4(0) = —4 a® (3 — 15 9 + 4 7) = a5 cos?6 (75 — 260 n + 96 77?) 
+ 2 a4 f’ cos@ (52 — 187 9 + 56 7?) + a4 f’ cos®6 (77 — 264 1 + 96 772) 
+ 8 a f cos (2 — 8 + 3 9%) — a3 f” (37 — 122 n + 32 7?) (3.6) 
= 16 a f” cos?6 (9 — 32 + 10 92) + a2 f’”* cos (87 — 314 7 + 88 7?) 
—2af” (8-—31y+87)), 


Q;(0) = (1 — n) a{2 (2+ 9 9 — 4 9? — 16 7) I(c — 2, s, 1) 
(4 — 7 9 = 30 7 + 8 4? + 32.y5) s- [le = 2, 5, 1)} 
+ (3 — y) a{— 2 (22 — 47 9 + 44 7? — 16 7) I(c, s, 1) 
4+—(32 + 11 » — 82°77 + 88 1° — 32-4) s-* Je, s, 1} 
+- f cos0 {2 (39 — 109 n + 162 y? — 120 3 + 32 7) I(c—1,s, 1) 
(27 — 27 y — 128 7? + 300 73 — 240 74 + 64 75) 
Kis tef (a= Ipseat 


(3.7) 
— (1 — n) a? f’ cos6 {2 (1 — 4 7) I(c — 1,s, 3) 
~ (5+ 24-89) s+ J(c— 1,5, 3)} 
+ (3 — ») a? f’ cos6 {2 (5 — 4) I(c +1, s, 3) 
+ (3 —107+8 7?) s1 J(c+1,s, 3)} 
+ 2 a8 I(c, s, 3) {— 44 + 61 y — 28 7? + (24 — 35 4 + 20 17?) cos?6} 
+ a®s-4 J(c, s, 3) {— 18 + 81 4 — 122 7? + 36 7° 
+ (6 — 35 4 + 70 7? — 40 77) cos?6} , 
and 
Q,(0) = a® (5 — 48 7 + 32 7?) — 8 a5 cos? (22 9 — 25 7? + 8 77°) 
— at f’ cosO (23 — 306 9 + 232 4° — 32 7°) 
+ 4 a f’ cos? (45 4 — 52 7? + 16 99) + 8 a4 f cos6 (5 4 — 6 77? + 2 97°) BA 


4 a3 f” (13 — 114 7 + 64 972) + 2.43 f” cos®O (13 — 200 + 172 4? — 32%) 
— a? f” cos (29 — 266 7 + 184 7? — 16 7) 
4+2af” (4—27y + 167%). 
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In equation (3.5) we have 


R, cos0, = acos6 — f’, 


and it may be noted that Q,(8) is finite on the axis 6 = 0. We define the integrals 
I(m, s,n) and J(m, s,) which occur in equations (3.4) and (3.7) as 


1 ili 
> xm logé “2” sin (s log 
I(m, s, n) = j= eee di, J(m,s,n) =| a 84) dd , 
J A 
0 0 


where Ri = x2 + y?+ (z—A/’)?. Further, we have 4s?=3—47? and 
c = 3/2 — 7. 

Since the expressions (3.1) and (3.2) do not lend themselves so readily to 
numerical computation as the expressions for the stresses in series of spherical 
harmonics, we now give these latter expressions. These can be found by using 
equations (2.2) to (2.5) together with equations (3.16), (3.18), (4.5) and (4.7) of 
SH. An alternative method is to obtain solutions of the equations satisfied by 
the displacement components as series of spherical harmonics suitable for the 
region r > a, the coefficients of the harmonics in these series being found from 
the condition that the boundary v = a of the hollow is stress free. The additional 


stresses 77, and r 6, are given by 


ii = — 


2 eos (m — 1) (n?+ 3n— 27) eae 


(af)? — (2m — 1) W(n) y 
(3.9) 
(1 +1) U,(n) (a\n+3] (ayn 
(2n + 3) W(m) () (7) F,, (cos) , 
and 
ie et — (a —1) (n?— 2+ 2n) Un) (a\nrt 
Lee hemi eagtes 
(3.10) 


U,(n) A\n+3 a\n 
~ ena ay wen Cr) | (FG) Zi ose), 
where 
W(n) =n? + (L1—2n)n+1—yn, 
Ug(n) =? (2% +1) (n+ 3 —4 7) —a*(2n—1) (n +1), 


U,(n) = m (n — 1) (w + 2) [f? ( 2n 4+ 3) (1 + 3-47) — a (2+ 1) (w+ 1)] 


—4a(2n+1) (1-7). 
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The stresses 6 6 and ¢@ on the boundary r = a are obtained as 


coos ae 
6 7) = Bac » Ne - Vo(n) : a\n ’ 
(a f)* ant (2 ” + 3) (n + 2) W(n) (+) Pe (cos 0) | 
co (3.23) 
2m A (1 — ») coté Vi (n) Pana 
ss (a f)* (2 + 3) (n + 2) W(n) i T; (cos6) , 
and 
Ry 2uaA ~ , V,(n) fa\n 
PP “af? Pm (2 + 3) (n+ 2) W(n) (+) P, (cos 6) | 
Sel 
_ 2A (1 — 9) cotd y ra) a\" or ears 
(a f)? = (2 + 3) (n + 2) W(n) (|) , (cos) , 


where 


Vo() = a? (m + 1) [4 05 + 4 (4 — y) nt + (19 — 16 ) 2 


r (14 —19 9 — 8 9?) n? + 2 (8 —3n — 872) n+ 6 (1 — 7) 
—f? (2+ 1) (2 +3)" (n+ 2) (n+3—4%) (n—7), 
V,(n) = a? [4 n* + 20 n3 + (27 — 8 y) n? + (9 — 16 7) xn — 64) 
— f? (2 + 1) (2”+ 3) (n+ 2) (w2+3—4y) 


Vi (n) = a® (n + 1) [44 n° — 4 (1 — 4) n4— (8-199, + 8 2) 23 
+ (5+ 6 — 24 7) n? 4+ 248 — 11 7) n + 6 (1 — 7)] 
— f? (2+ 1) (2n+ 3) n (n+ 2) (n+3—4y) (yn—1). 


To illustrate the variation of the stresses 6 6 and @ @ on the boundary r = a 
these stresses have been calculated for different values of f and 7. These 
calculations were performed by truncating each of the series in (3.11) and 
(3.12) after N terms, N being chosen so that the truncation error in each case 
is sufficiently small to ensure the calculated values of the stresses are correct to 
four decimal places. The sums of the series to N terms were then calculated on 
the University of Durham Computing Laboratory Pegasus Computer using a 
method due to CLENSHAW [8]. The values of the stresses thus obtained are 
given in Tables 1 and 2 and plotted in Figures 1 to 3. The actual values shown 


are those of 6 6lu A and G@/u A. 
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Table 1 


ZAMP 


Elastic Solid Containing a Spherical Hollow. Values of 6 Olu Aonv=4a 
Oe eS SS SS 


n = 0-2 H = O-4 

7) 

| jee f=4 = 10 j=2 

0° — 9-2069 — 0-8371 — 0:0814 — 9-9313 
10° 5-7950 0-6857 0-0709 6:6079 
20° +0-5016 0:2950 0-0417 0-4508 
307 4-6533 + 0-1871 + 0-0003 + 3:6740 
40° 5-9566 0-6148 0:0472 5:0752 
50° 5-6344 0-9030 90-0909 4-:9113 
60° 4-7072 1-0334 0-1252 4-1478 
70° 3-6822 1-0308 0-1461 3-2651 
80° 2°7550 0:9357 0-1526 2:4509 
90° 41-9795 0-7870 0:1458 1-7604 
100° 1-3545 0-6156 0-1280 1-1971 
110° 0-8611 0:4427 0-1026 0°7471 
120° 0-4771 0-2818 0-0731 0-3930 
130° 0-1828 0-1407 0-0426 0-1188 
140° — 0-0380 0-0234 0-0139 — 0-0889 
150° 0-1979 — 0-0683 — 0-0106 0-2405 
160° 0-3060 0-1337 0-0293 0-3436 
170° 0-3685 0-1729 0-0410 0-4034 
180° 0-3889 0-1859 0-0449 0-4231 


Figure 1 


Elastic solid containing a spherical hole. 
Graphs of 66/4 A and ¢@/mu A for f = 2, 
n = 0-2. — 96] A; —--Golm A. 


Figure 2 
Elastic solid containing a spherical hole. 
Graphs of 6 O/u A and ¢@/u A for f = 2, 
4 = 0-4. —66/u A; ---GQlp A. 
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Table 2 
Elastic Solid Containing a Spherical Hollow. Values of PG PluA onr=a 


- H = 0:2 | y = 0-4 

| f=2 eee f—10 Pee 

OP — 9-2069 —0-8371 — 0-0814 —9+9313 
10° 7:9420 0-7940 0-0790 8-1241 
20° 5+3457 0-6801 0-0721 4.5680 
30° SATUS 0-5318 0-0620 1-7761 
40° 17288 0-3866 0-0505 0-3088 
50° 0-9923 0-2685 0-0393 + 0-2694- 
60° 0-6275 0-1861 0-0298 0-4042 
70° 0-4560 0-1367 0-0230 0-3513 
80° 0-3811 0-1131 0-0190 0:2327 
90° 0-3533 0-1074 0-0179 0-1019 
100° 0-3477 0-1125 0-0191 —0:0196 
110° 0-3518 0-1236 0-0221 0-1245 
120° 0-3595 0-1370 0:0263 0-2112 
130° 0-3678 0-1504 0-0308 0-2806 
140° 0-3752 0-1626 0-0353 0-3344 
150° 0-3812 0-1726 0-0393 0-3743 
160° 0-3855 0-1799 0-0423 0-4018 
170° 0-3881 0-1844 0-0443 0-4178 
180° 0-3889 0-1859 0-0449 0-4231 


It is of interest to compare the variation of the stress 9 6 in this problem 
with the variation of the hoop stress round the boundary of a circular hole in 
an infinite plate in a state of generalized plane stress due to the action of a 
force nucleus at a point of the plate, the nucleus being directed radially out- 
wards from the centre of the hole. This problem has been discussed by HOLGATE 
[9] and the stress curves given by him are found to be similar to the curves 
given in Figures 1 to 3. 


Figure 3 
Elastic solid containing a spherical hole. Graphs of 66 Olu A and 1p Plu A for f = 4 and f = 10, 
y = 0-2 (note different scale from Figures 1 and 2). —0 A Olu A} —==—©@ ole A, 
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4. Force Nuclei Acting at Internal Points of an Elastic Sphere 


We now consider the problem of an elastic sphere under the action of two 
equal and opposite force nuclei, magnitude F, applied at points on a diameter 
of the sphere equidistant from the centre, the axes of the nuclei being directed 
along this diameter. Since the case when the forces act at the surface of the 
sphere has been discussed by STERNBERG and ROSENTHAL [10], we suppose the 
nuclei to act at internal points, the boundary of the sphere being unstressed. As 
the method given in SH for an infinite solid containing a spherical hollow can 
be extended to an elastic sphere (see SH, p. 250), it is possible to find closed 
form expressions for the stresses similar to those found in section 3. Since it 
is again more convenient for numerical work to use expressions for the stresses 
as infinite series of spherical harmonics, we give only these expressions. 

If the nuclei act at the points whose cartesian coordinates are (0, 0, /) and 
(0,0, —f), (f < a), and are directed radially outwards from the centre of the 
sphere, the stresses due to these nuclei acting in an infinite unbounded elastic 
solid are, when f < 7, 


7}, = —*08 eee (4n?> + 6n—2 )) eae 


n=0 (4 — 1) Y 


(2+ 1) (2+ 2) (2n—244n) / f\2n+8 
+e ee eae 


705 = ey ee 2 n) apes 


Pe = (4 — 1) F: 
(2m + 2)(2n—-2+4n) /f (4.2) 
Zn — 2n+3 
Sree era eh pe 
Af _ 4HA * 2n(4n®?—4n—1-+4 27) f \2n+2 
06, Vie Pa) (4n— 1) ] (=) 
_ (2n4+ 1)? (2n—2+47n) / fants 
feta le) | Pamteost) 
(4.3) 


4 4 uA cot Se [2 7a Baa 


(2m — 2+ 47) f\2n+3 
(4 + 3) (5) ‘| T;,, (cos 6) 


(4 
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and 


_. (2m+1) (2n—24 470) (fy"""] P,,, (cos) 


(4 n + 3) Y 4 4) 
4 4uA cot 6 (2n—4+ 4%) /( f\2n+1 | : 
‘ia ae [- 4 u— ibs =) 
(2m—2+4n) / f\2"+3 | 
<a Cy) | Ten 0080), 


where F = 16” uA (1—7). When r < f, the expressions for the stresses are 
obtained from equations (4.1) to (4.4) by replacing » by — (nm + 1) in the n-th 
term of each series. 

The additional stresses 77, and 70, which must be added to the stresses 77, 


and 70, respectively to give the stress distribution due to the force nuclei 
acting in the elastic sphere, can be found by methods similar to those given 
in section 3 for the solid containing a spherical hollow. We find that 


rh, = oes silent od. SF ede | 


(af? — (4.n + 3) C(n) a ae 
n A,(n) vy \2n-2 vi Q2n+1 | 
+ (an — 1) Ci) (Sy | GG) Bee (0088) 
and 
ee 8uA oe (n+ 1) (4n?+4n—1+ 2m) A,(n) (7\2" | 
r= apr & | (tn + 3) C(n) i) es 
A,(n) y \2n—2 ii 2n+1 4 | 
(4m — 1) C(n) (5) | (1) Toy (COs) , 
where 


Cin) =4r74+2(14+2n)n+1+%, 
A,(n) = a2 (4n+3)n—f (4n+1) (n-14+27), 


A,(n) = (n + 1) (4? — 1) [a2 (4n+1)n—f? (4n—1) (n—-1+27y)] 
+ a2(4n+1) (1-7). 
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Further, we obtain the stresses 06 and @@ on the surface of the sphere as | 


00= me yy Deve tan 1) C(n) ake Fan (ost) r 
8 cee Ie . SS om wae TCG € J" Fay (C088), 

and 

= a oy (2n— sien 1) C(n) ae "Pay (0080) | d 
Se 2 CLES) (7) """ Ti (6086) | 


B,(n) = a? n [128 n5 + 64 (1 + 7) nt — 40 n3 — 4 (13 — 8) (1+ 9) 
+2n+ (3-27) (+n) 
— f2 (16 n? — 1) (4n2—1) (2n 4147) (w—-14+27), 
B,(n) = a? [32 n* — 16 "8 —2(9+ 87) w?+n+1 +4 7) 
— f? (16 m2? — 1) (2n—1) w»—142n), 
and 
B,(n) = a? n [128 n n® + 64 (1 + 7) n4*+ 8 (8 —5 y — 8 7?) n8 
— 20 (14 9) ?-—2(2—n-—27r?)n+1+7) 
= 7? (16 0? = 1) (4.8 — D295 ed ie = ee 


To illustrate the variation of the stresses 6 6 and @ over the surface of the 
sphere, these stresses have been calculated for different values of f and y using 
equations (4.7) and (4.8). The values obtained are given in Tables 3 and 4 
and plotted in Figures 4 to 6, the actual values shown being those of 0 6/8 A 
and ¢ g/8 uw A. Since the stresses have mirror symmetry about the plane z = 0 
only the values for 9 between 0° and 90° are given. 


In conclusion, I wish to thank the University of Durham Computing 
Laboratory for the facilities placed at my disposal. 


, 
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Table 3 
Elastic Sphere. Values of 66/8 eAonr=a 


= 0-2 n= Or’ 


0 ———* 
| el | f= 0:3 | 1208 | 40K 
0° + 0:5834 + 2-1196 + 5-4008 + 5-0576 
10° 0-5582 1-9592 44333 4-1131 
20° 0-4867 1-5469 2:5456 22776 
30° 0-3796 1-0345 1-0773 0-8630 
40° 0-2523 0-5580 0:2966 0-1241 
50° 0-1214 0-1858 —0-0421 —0-1855 
60° 0-0027 —0-0717 0-1695 0-2939 
70° —0-0911 0-2322 02103 0-3230 
go° 0-1509 0-3178 0-2203 0-3266 
| 90° 0-1715 0-3445 | 0-2217 0-3259 
Table 4 
Elastic Sphere. Values of Gp|/8 uA onr =a 
9 y = 0:2 nH = 0-4 
f= O41 = 0:3 f = 0-5 f= 0-5 
0° + 0-5834 + 2-1196 +5-4008 +5-0576 
10° 0-5580 1-9833 4-7091 4.2918 
20° 0-4855 1-6205 3-2059 2-6901 
30° 0-3761 1-1384 1-7661 1-2588 
40° 0-2445 0-6464 0-7366 0-3223 
50° 0-1076 0-2170 0-0859 —0:2151 
60° —0-0182 —0-1170 — 0-3035 0-5075 
70° 0-1188 0-3493 0-5259 0-6607 
80° 0-1835 0-4846 0-6396 0-7340 
90° 0-2058 0-5288 0-6745 0-7557 


4 Figure 4 
Elastic sphere. Graphs of 6/8 uA and @ @/8 A for f = 0-1 and f = 0:3, 7 = 0-2. 
— 66/8" A; ---G@/8 MA. 
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Qo 


Figure 6 
Elastic sphere. Graphs of 00/8 uA and 
PQ/8 “A for f = 0-5, y = 0-4 : 
(note different scale from Figure 4). 


Figure 5 
Elastic sphere. Graphs of 0 0/8“ A and 
PQS “A for f = 0-5, y = 02 
(note different scale from Figure 4). 


REFERENCES 


(1) R. D. Minprin, Physics 7, 195-202 (1936). 
[2] B. Sen, Bull. Calcutta math. Soc. 22, 73-83 (1940). 
[3] W. R. Dean, H. W. Parsons and I. N. SNEDpon, Proc. Camb. phil. Soc. 40, _ 
5-19 (1944). 
[4] W. D. Cottins, Mathematika 7, 125-130 (1954). 
[5] W. D. Cotttns, Proc. London math. Soc. 9 (3), 242-252 (1959). 
[6] A. E. H. Love, Mathematical Theory of Elasticity (Cambridge, 1927), p. 185-186. — 
[7] T. M. MacRopert, Spherical Harmonics (Methuen, 1947), p. 125. 
[8] C. W. CLensHaw, M. T. A. C. 9, 118-120 (1955). 
[9] S. Horeate, Proc. Camb. phil. Soc. 40, 172-188 (1944). 
[10] E. STERNBERG and F. ROSENTHAL, J. appl Mech. 79, 413-421 (1952). 


Résumé 


Dans cette communication, nous donnons des expressions pour les efforts causés 
par des noyaux de forces agissant a certains points 

(i) d’un solide élastique contenant une cavité sphérique; 

(ii) d’une boule élastique. 
Dans le premier cas, on obtient des expressions en séries de fonctions de LEGENDRE. 
Des expressions équivalentes contenant des intégrales et des fonctions élémentaires 
sont également obtenues. Dans le deuxiéme cas, on obtient seulement des expres- 
sions en fonction de LEGENDRE. Dans chaque cas, des résultats numériques sont 
obtenus pour les efforts sur la surface sphérique. 


(Received: July 17, 1959.) 
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Uber die Erzeugung von hohen, kurzzcitigen Magnetfeldern 


Von Prero Corti, Ziirich?) 


1. Einleitung 


Die Erzeugung von Dauermagnetfeldern tiber 60000 A/cm (75000 Oe) er- 
fordert bekanntlich einen so grossen apparativen Aufwand, dass nur sehr wenige 
Laboratorien iiber die entsprechenden Installationen verfiigen. Die héchsten 
Dauermagnetfelder, die bis jetzt hergestellt werden konnten, liegen bei etwa 
80000 A/cm (100000 Oe). Wie Kapitza [1]?) aber schon 1924 gezeigt hat, ist 
es mit einem relativ bescheidenen Aufwand méglich, drei- bis viermal so starke 
Magnetfelder zu erzeugen, sobald man sich darauf beschrankt, diese nur fiir 
ein sehr kurzes Zeitintervall (Gréssenordnung Millisekunden) zur Verfiigung zu 
halten. 

KapiItzA verwendete zuerst eine Serie von Akkumulatoren als Energie- 
quelle. Mit einem mechanischen Schalter wurde diese wahrend einiger Hundert- 
stelsekunden iiber eine kleine Spule kurzgeschlossen, wobei Magnetfelder bis 
zu 100000 A/cm erzeugt wurden. Um noch héhere Magnetfelder zu erreichen 
und diese auch in einem grésseren Volumen von mehreren Kubikzentimetern 
zur Verfiigung zu haben, beniitzte KAPITzA [2] spater seine beriihmt gewordene 
Dynamomaschine als Energiequelle. Er konnte so wahrend einer Impulsdauer 
von rund 30 ms in einem Volumen von ungefahr 6 cm? Magnetfelder bis zu 
256000 A/cm (320000 Oe) erzeugen. 

Eine starke Reduktion des apparativen Aufwandes brachte 1926 die Idee 
von WALL [3], eine Batterie von Kondensatoren als Energiequelle zu beniitzen. 
Die Kondensatorenbatterie von WALL wies eine Kapazitat von 1400 pF auf 
und konnte auf maximal 2000 V aufgeladen werden. Durch die Entladung die- 
ser Kondensatorenbatterie iiber eine kleine Hochfeldspule konnte WALL Ma- 
gnetfelder bis zu 182000 A/cm (228000 Oe) erzeugen, wobei allerdings die 
Hochfeldspule infolge der Biot-Savart-Krafte explodierte. Diese Technik wurde 
seither in fast allen Laboratorien [4-10], die sich mit der Erzeugung von kurz- 
zeitigen, hohen Magnetfeldern beschaftigt haben, verwendet. 

Im Verlaufe der letzten 5 Jahre sind am Institut fiir kalorische Apparate 
und Kiltetechnik der ETH einige Apparaturen gebaut worden, die ebenfalls 
nach diesem Prinzip arbeiten. Diese Apparaturen wurden hauptsdchlich fiir die 
Messung der Anderung des elektrischen Widerstandes von reinen Metallen in 


1) ETH. Institut fiir kalorische Apparate und Kaltetechnik. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 32. 
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hohen Magnetfeldern beniitzt. Dazu sind relativ lange Impulszeiten erforderlich, 
denn zu schnell veranderliche Felder werden im allgemeinen die gut leitenden 
Metallproben nicht geniigend durchdringen und sie auch durch starke Wirbel- 
stréme aufheizen. Es wurden hier zum Beispiel wahrend einer Impulsdauer von 
ungefahr einer Millisekunde in einem Volumen von 0,5 cm’ Magnetfelder bis zu 
320000 A/cm (400000 Oe), oder wahrend etwa 100 Millisekunden und in einem 
Volumen von 10 cm Magnetfelder bis zu 80000 A/cm (100000 Oe) erzeugt. 

Uber die Erzeugung noch héherer Felder (bis zu 750000 Oe) wahrend einer 
Zeit von gréssenordnungsmassig 100 ys sind in den letzten Jahren drei ausfiihr- 
liche Arbeiten erschienen [6, 7, 18]. Erganzend soll hier im folgenden die Er- 
zeugung von Hochfeldern wahrend relativ langer Zeiten und die 6konomisch 
giinstigste Dimensionierung einer solchen Hochfeldapparatur behandelt werden. 
Die hier angefiihrten Berechnungen folgen zum Teil einer wertvollen Arbeit 
von CHAMPION [4], sind aber speziell unseren Zwecken angepasst worden. 

Die Messmethode, die hier angewendet wird, soll am Schluss kurz beschrie- 
ben und eine Messung der magnetischen Widerstandsanderung von Kupfer ge- 
zeigt werden. 


2. Theorie 
Figur 1 zeigt das prinzipielle Schaltbild unserer Apparatur. Die Strom- 
starke in einem solchen Entladungskreis wird 


U =Z D) . 
I(t) — eee é URS sIn@w t ) (1) 


raed Fide ee 
ra 0 ahr (2) 
U bezeichnet hier die Aufladespannung der Kondensatorenbatterie, L die 


Selbstinduktion der Hochfeldspule, C die Kapazitat der Kondensatorenbatterie 
und R den Ohmschen Widerstand des Entladungskreises. 


wobei 


Schalter-Steuerung 


Aufiade - 
Gerat 


Figur 1 
Prinzipielles Schaltschema der Hochfeldapparatur. 
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Diese Beziehung gilt sowohl fiir reelle als auch fiir imaginare w und geht im 
aperiodischen Grenzfall wm > 0 gegen den Ausdruck 


U Fe 
T(t) be SBRht (3) 


Wie die Erfahrung zeigt, liegt die Art unserer Entladungen meistens in der 
Nahe dieses aperiodischen Grenzfalles. 


C U E H V 
(pA IV) [J] [kOe] [Alem] — [cm] 
iT 5 +198 
105 10 
0 = ee = yg? 
4-797 ‘iN 
she 10° 10 ate 
107+ 10° eG er a ant 
+f sla ea eee / 10 
Bree 3 Cae es 4 
wt -- oe : ® ss! 
Sy 4 210° si 
107 10 zx 
ott 10? +19 
iL : ii 
‘ 10 10 10 
10 10 fa’ 
; 10° 


Figur 2 
Nomogramm. Das maximal erzeugbare Magnetfeld H als Funktion der eingesetzten Energie E und 
des inneren Volumens V der Hochfeldspule (energetischer Idealfall). Zur Auswertung sind die drei 
linken und die drei rechten Funktionsskalen je durch eine Gerade zu verbinden (C und U sind 
Kapazitat und Aufladespannung der Kondensatorenbatterie). 


Im energetisch giinstigsten Fall wiirde die gesamte, urspriinglich in den 
Kondensatoren aufgespeicherte, elektrische Energie E = C U?/2 im Augenblick 
des Feldmaximums lediglich in Form von magnetischer Energie im innern 
Volumen der Hochfeldspule stecken. In Giorgi-Einheiten wiirde also 

CU2 MA; 


E=—7- =“ 7, (4) 


wobei H,,,, das Feldmaximum im Mittelpunkt der Spule und V das innere 
Volumen der Spule bezeichnet. 

Das Nomogramm von Figur 2 gibt eine Darstellung dieses Zusammenhanges. 
_ Wir werden allerdings sehen, dass der Wirkungsgrad dieser Energieumformung 
selten iiber 50% liegt. Es ist zweckmassig, diesen Wirkungsgrad mit K® zu be- 
zeichnen. Dann wird 


Ena = KU YS, Que eee ay (5) 


Haz [A/em], C [uF], V [cm], U [V}. 
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Nach CHAMPION kann K in zwei Faktoren aufgespaltet werden: 
o<J<i,| 


(6) 
(ese r| 


R ekls 


Der Faktor J beriicksichtigt die Ohmschen Energieverluste, wahrend S der 
Tatsache Rechnung tragt, dass in einer Spule die magnetische Feldenergie nicht 
nur im Spuleninnern konzentriert werden kann. CHAMPION, der sich ausfiihr- 
lich mit diesen beiden Faktoren auseinandergesetzt hat, gibt einen expliziten 
Ausdruck fiir S. Tabelle 1 wurde einer seiner graphischen Darstellungen ent- 
nommen. 


Tabelle 1 
Der Faktor S als Funktion von « und B (siehe auch Figur 4) 


- ileal 1,5 2 4 8 
sce 
0,5 0,60 0,50 0,40 0,23 | 0,12 
1 0,82 0,70 0,57 0,32 0,17 
2 0,96 0,83 0,70 0,41 0,21 
4 0,98 0,88 0,75 0,47 0,26 
8 0,99 0,88 0,75 0,48 0,28 
Wie man sich leicht iiberlegt, gilt fiir den Faktor J 
1h 
Porat (7) 
x* (R=0) 
Mit der Substitution 
eth 
Yaar: (8) 
erhalt man aus Gleichung (1) 
if =s e~ (lVy=1) arctg Vy—1 (9) 


: ee monoton mit y (Figur 3). y = 1 entspricht dem aperiodischen Grenz- 
aul, 

Man wird also eine Hochfeldspule so dimensionieren, dass y méglichst gross 
und zugleich S nicht zu klein wird. Im Bereiche kleiner y-Werte ist vor allem 
die erste Bedingung massgebend, wahrend fiir y-Werte iiber 5 mehr Gewicht 
auf die zweite Bedingung gelegt werden muss. Ohne eine Nebenbedingung zu 
stellen, wird man dabei auf den Extremfall einer Spule gefiihrt, die aus einer 
einzelnen Windung von méglichst grossen Dimensionen besteht. 
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Jae Vp eeavet 


; y 
= ay = ! l i J bee 1 1 eb Set 
10 10 1 10 10? 10° 


Figur 3 


Der Faktor J als Funktion von y. 


Als Nebenbedingung soll die Selbstinduktion L der Hochfeldspule vor- 
geschrieben werden, da L im wesentlichen die Impulsdauer der Entladung be- 
stimmt. Fiir eine Spule von rechteckigem Wicklungsquerschnitt gilt bekannt- 
lich (Bezeichnungen nach Figur 4) 


Lon 7, Xa, 8), (10) 


wobei L in Mikrohenry und 7; in Zentimeter einzusetzen sind, wahrend / einen 
Formfaktor bezeichnet, der nur von den relativen Spulenparametern « und £ 
abhangt (a = 7,/7;, B =1/r;, siehe Tabelle 2). 


Tabelle 2 
Der Formfaktor 1, L = n? 7,1, L in Mikvohenry, vr; in Zentimeter 


| 
ax 
1,1 1,5 2 4 8 
tie ogee | 
0,5 24,5 23,0 25,5 B55 65 
3 14,0 16,5 19,0 SGD 58 
2 8,5 10,3 L239: 24,0 47,8 
4 4,8 6,2 7,6 15,8 37,0 
8 2,55 3,2 4,0 9,3 24,0 | 
Der Ohmsche Widerstand R einer Spule wird 
2 
qk on? x (a + 1) (11) 


~ 2fre—1) 8’ 


-wobei / den Fiillungsfaktor nach Figur 4 und @ den spezifischen Widerstand 
bezeichnet. Die Tatsache, dass sich 9 wahrend einer Entladung infolge der Er- 
warmung der Hochfeldspule verandert, soll vorlaufig vernachlassigt werden. 
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Wenn die Windungszahl » durch die Selbstinduktion L substituiert wird, er- 
halt man fiir y 


162) % 
ae ee ee oe 
wobel 1 2 
F(a, B) = (<= ; BA) (13) 


Figur 4 
Bezeichnungen der Hochfeldspule: r,/r; = «%; l/r; = B;  Windungszahl; @ Drahtdurchmesser; 
f Fiillungsfaktor = n 2227/8 1 Gy) e 


Der Geometriefaktor F wachst sehr schnell mit « und etwa 1000mal langsamer 
mit f (Tabelle 3). 


Tabelle 3 
Der Geometriefaktor F(a, B) in 10-18 Hy?/cm? 


Die energetisch giinstigste Dimensionierung der Hochfeldapparatur kann 
nun aus Gleichung (12) herausgelesen werden. Die Kondensatorenbatterie soll 
so angelegt werden, dass sie bei relativ kleiner Kapazitat auf mdglichst hohe 
Spannungen aufgeladen werden kann (Serieschaltung). Die héchsten Auflade- 
spannungen werden durch den Entladeschalter, die Isolation der Hochfeld- 
spule und nicht zuletzt auch durch Sicherheitsbedingungen beschrankt. Da 
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ausser der Kapazitat C nach (12) auch die Selbstinduktion L moglichst klein 
sein sollte, ist es grundsatzlich méglich, bei kiirzeren Impulsdauern einen gros- 
seren energetischen Wirkungsgrad zu realisieren. Z darf aber auch aus energeti- 
schen Griinden nicht zu klein gemacht werden, denn die Selbstinduktion der 
Hochfeldspule muss gréssenordnungsmissig grésser als diejenige der Zulei- 
tungen und der Kondensatorenbatterie sein, da sich ja die magnetische Energie 
im Verhaltnis der Selbstinduktionen aufteilt. 

Gleichung (12) besagt weiter, dass eine beliebige Verkleinerung der Hoch- 
feldspule (v;) bei gegebener Kondensatorenbatterie nicht unbedingt ein hdhe- 
res Magnetfeld liefert, da y mit der vierten Potenz von 7; zunimmt und J im 
Bereich kleiner y-Werte ebenfalls sehr rasch mit yy wachst. 

Der Fiillungsfaktor f tritt in Gleichung (12) im Quadrat auf. Da sich ein 
grosses /, wie man sich leicht iiberlegt, auch giinstig auf die mechanische Wi- 
derstandsfahigkeit einer Hochfeldspule auswirkt, ist es in doppelter Hinsicht 
vorteilhaft, Rechteckprofildraht zu verwenden. 

Im allgemeinen wird das innere Volumen der Hochfeldspule und die mini- 
male Impulsdauer der Entladung durch die auszufiihrenden Messungen vorge- 
schrieben. Anhand der Thomsonschen Formel kann aus der minimalen Impuls- 
dauer ein Wert fiir L festgelegt werden: 

4 Tih 

L= aeeea (14) 
wobei T,,,, den Zeitpunkt des Feldmaximums bezeichnet. Bekanntlich ist die 
Periodendauer von gedampften Schwingungen immer etwas grésser als der 
nach der Thomsonschen Formel berechnete Wert. Wie wir gesehen haben, 
nimmt y mit « schnell zu. Eine Vergrésserung von « hat allerdings eine Ver- 
schlechterung der Homogenitat des Magnetfeldes langs der Spulenachse zur 
Folge. Hat man einmal das grésste «, das sich diesbeziiglich noch verantworten 
lasst, bestimmt, so kann daraus anhand von (10) die Windungszahl n berechnet 
werden. Bei der Festlegung von « muss ausserdem darauf geachtet werden, 
dass der Faktor S mit zunehmendem « abnimmt. Eine Vergrésserung von « ist 
deshalb vor allem dann zweckmassig, wenn der y-Wert noch kleiner als 5 ist. 


3. Die Erwarmung der Hochfeldspule 


Da die Erwarmung einer Hochfeldspule im wesentlichen adiabatisch ver- 
lauft, gilt hier 
Pod = 0, Pal , (15) 


wobei 7 die Stromdichte, c, die spezifische Warme pro Gramm und » die Dichte 
bezeichnet. 
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Wir trennen die Temperatur- und die zeitabhangigen Grossen und integrieren 


t Tt 

ee Cy? 
pe a= | “oat, (16) 
0 T4 


wobei 7, die Anfangstemperatur und 7; die Temperatur zur Zeit t bezeichnet. 

Fiir Temperaturen unter 200° K kann das Integral rechts nur numerisch 
ausgewertet werden. OLSEN [5], der seine Hochfeldspulen in fliissiges Helium 
getaucht hatte, gibt eine graphische Darstellung dieses Zusammenhanges fiir 
Kupferspulen von 1,6 - 10-§ Q cm bzw. 0,8 - 10-8 Q cm spezifischen Restwider- 
stand. Diese graphische Darstellung geht bis zu Endtemperaturen von 100° K. 
Bei héheren Temperaturen ist c, ndherungsweise konstant und @ proportional 


zur absoluten Temperatur. Dann wird 
t 


Q273 . 
T; aegret 
fine l a7 
Anhand dieser Formel wurden in Figur 5 die graphischen Darstellungen von 
OLSEN erweitert. 
°K if 
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Ss 
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Figur 5 


Die Temperaturerhéhung der Hochfeldspule als Funktion von fe dt fiir drei verschiedene Aus- 
gangstemperaturen, T 4. 


Das Integral fi dt kann durch das Magnetfeld H ausgedriickt werden, indem 
fiir H der Naherungswert In/2 1 eingesetzt wird. Man erhalt so 


(i 12 ore He 2 li sys p T max 


(Yq os ¥,)? : 


wobei pT,,,, ein Mass fiir die Dauer der Entladung darstellt. Fiir eine ideal 
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sinusférmige Entladung (R = 0) wird p = 1, wahrend im aperiodischen Grenz- 
fall = 1,85 wird, 

Eine kleine Hochfeldspule (Tabelle 4, Nr. 1) wiirde sich zum Beispiel bei 
der Erzeugung von 250000 A/cm wahrend 3 ms von 80° K auf etwa 300° K, von 
300° K aber auf rund 1000° K erhitzen. Fiir die Erzeugung von Magnetfeldern 
iiber 200000 A/cm (wahrend Millisekunden und mehr) ist deshalb eine Vorkith- 
lung der Hochfeldspule mit fliissiger Luft unumgiénglich. Aus dem letzten 
Ausdruck geht jedoch hervor, dass sich die Temperaturerhéhung durch eine 


HSH 


Figur 6 


Die Messung der Temperaturerhéhung der Hochfeldspule. Die untere Kurve reprisentiert den 

Feldverlauf (100 Skt = 2440 A = 468000 A/em = 608000 Oe), wahrend die obere Kurve den 

Verlauf der Ohmschen Spannung an der Hochfeldspule darstellt (100 Skt = 1270 V). Der ganzen 
Bildbreite entspricht eine Zeit von 10 ms. 


Vergrésserung der Wicklungsstaérke wirksam reduzieren lasst. Um die Aus- 
driicke iiber die Temperaturerhéhung zu kontrollieren, wurde mit einem Ka- 
thodenstrahloszillographen der Verlauf des Ohmschen Widerstandes einer in 
fliissiger Luft vorgekiihlten Spule wahrend einer Entladung gemessen. 

Dabei wurde die induktive Komponente der an der Spule auftretenden’ 
Spannung mit einer Pick-up-Spule auskompensiert. Die Verschiebung der 
Maxima demonstriert augenfallig die Vergrésserung des Widerstandes wahrend 
der Entladung. Nach 2,5 ms betragt der Ohmsche Widerstand 1,65 Q, was 
bereits einer Temperatur von 215°K entspricht. Der Fiillungsfaktor dieser 
Spule war 0,55 und die Wicklungsstarke 0,6 cm. Um anhand von Formel (17) 
auf eine Endtemperatur von 215° K zu kommen, muss hier fiir if 7? dt ein Wert 
von 5,8- 108 A2s cm-4 eingesetzt werden. Dies entspricht einem #T,,,, von 
2,1 ms (H,,,+ = 160000 A/cm). 


4. Die Verstarkung der Hochfeldspule 


Eine unverstarkte Kupferspule von rechteckigem Wicklungsquerschnitt 
wird auch bei giinstigster Dimensionierung bei Magnetfeldern zwischen 200000 
und 250000 A/cm durch die Biot-Savart-Krafte in radialer Richtung aus- 
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einandergerissen. Bei Impulsdauern in der Gréssenordnung von Millisekunden 
und weniger setzt also nicht die Erwarmung der Hochfeldspule der Erzeugung 
von Magnetfeldern zuerst eine Grenze! ee 
Fiir die azimutale Zugbeanspruchung o einer Hochfeldspule durch die Biot- 
Savart-Krafte gilt, unter Vernachlassigung der zusatzlichen Belastung der 
4usseren Windungen durch den Auflagedruck der inneren Windungen, 


2 Ma dA, LaPenites (19) 


wobei y den Abstand von der Spulenachse bezeichnet. Bei der Herleitung ~ 
dieser Formel wurde fiir H der Naherungswert I /2 / eingesetzt, der bekannt- — 
lich um so genauer stimmt, je langer die betrachtete Spule gegentiber ihrem 
mittleren Durchmesser ist. Die maximale Zugbeanspruchung tritt fiir Spulen, 
deren « grésser als 2 ist, bei y= 7,/2, und fiir Spulen, deren « kleiner oder 
gleich 2 ist (wo sie bei gleichem H grésser wird), beiy=7; auf. Fir ein H von zum 
Beispiel 200000 A/cm («= 2) erhalt man fiir die maximale Zugbeanspruchung ~ 
ungefahr 50000 N/cm?2, was ungefahr dem doppelten Wert der maximal zulas- — 
sigen Zugbeanspruchung von handelsiiblichem Kupfer entspricht. Die starkere — 
Beanspruchung der inneren Windungen einer Hochfeldspule hat nach meh- — 
reren Entladungen eine Ermiidung der Widerstandsfahigkeit zur Folge, indem 
bei jeder Entladung die inneren Windungen etwas mehr gegen die dusseren — 
Windungen fliessen, bis die Deformation im Innern zu gross und die Spule 
unterbrochen wird. Einlagige Kupferspulen kénnen schon bei 100000 A/cm 
explodieren. Eine unverstarkte Kupferspule, die Felder bis zu 200000 A/cm 
aushalten soll, muss deshalb eine méglichst grosse Wicklungsstarke (« gross) 
aufweisen und wenn méglich, zur Vergrésserung des Fiillungsfaktors f/, aus 
Rechteckprofildraht gewickelt werden. 

Von den verschiedenen Verstarkungsméglichkeiten, die hier ausprobiert 
wurden, haben sich verleimte Glasfasern am besten bewahrt’). Solche Glas- 
fasern sind in Glasseidenbandform im Handel erhaltlich. Wir haben ein 10 mm 
breites, etwa 0,05 mm dickes Glasseidenband verwendet, das aus den diinnsten 
serienmassig hergestellten Mikrofaden gesponnen wird, da so diinne Glasfasern 
die grésste spezifische Reissfestigkeit aufweisen. Dieses Band wurde in einem 
Giessharz (zum Beispiel Araldit) getrankt und dann moglichst straff in ein- 
zelnen Schichten zwischen die Wicklungslagen und um die Hochfeldspule 
gewickelt. 

Als Beispiel fiir eine solche Glasfaserverstarkung mége der Aufbau von 
Spule Nr. 1 (Tabelle 4) dienen. Diese Spule, die speziell fiir die Erzeugung 
hochster Magnetfelder wahrend einer Impulsdauer von rund einer Millisekunde 
dimensioniert ist, wurde aus lackisoliertem Rechteckprofildraht (1,30 - 0,25mm) 


*) Fiir diesen Vorschlag bin ich Herrn Prof. Dr. P. GRASSMANN zu Dank verpflichtet. 
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auf einen Plexiglaskérper gewickelt. Nach sechs Lagen wurde eine ungefahr 
15mm dicke Glasfaserschicht eingeschaltet. Dabei musste speziell darauf 
geachtet werden, dass sich die Glasfasern vollstandig mit Giessharz durch- 
setzten (das unverleimte Glasfaserband kann mit Leichtigkeit mit der Hand 
zerrissen werden). Nach 10 weiteren Lagen (total 231 Windungen) wurde noch 
eine etwa 5 mm dicke Glasfaserschicht um die ganze Spule gelegt. Nach einer 
Vorkihlung mit fliissiger Luft konnte in dieser Spule ein maximales Magnet- 
feld von 290000 A/em (360000 Oe) erzeugt werden, wihrend sie bei 
312000 A/cm (390000 Oe) explodierte. Es konnte allerdings nicht eindeutig 
festgestellt werden, ob die Zerstérung infolge der Biot-Savart-Krafte oder 
infolge der Erhitzung erfolgte. Der aussere Glasfasermantel ist jedenfalls dabei 
unversehrt geblieben. 

Durch den Ubergang zu kiirzeren Impulszeiten kénnten wahrscheinlich in 
solchen mit Glasfasern verstarkten Spulen noch betrachtlich héhere Magnet- 
felder erzeugt werden, denn einerseits kann dadurch die Temperaturerhéhung 
reduziert werden und andererseits spielt besonders bei sehr kurzen Impuls- 
dauern (zum Beispiel 10 us) auch die Tragheit des Spulenmaterials eine ent- 
scheidende Rolle, indem sich die Hochfeldspule wahrend einer so kurzen Zeit 
gar nicht geniigend deformieren kann, um zu zerreissen (siehe zum Beispiel [7]). 


5. Die Hochfeldapparatur 
a) Kondensatorenbatterie 


Unsere Kondensatorenbatterie besteht zurzeit aus 80 Elektrolytkonden- 
satoren yon je 550uF Kapazitaét und 500 V Nennspannung (total 5500 J 
Fassungsvermégen). Es wurden vor allem aus finanziellen Griinden Elektrolyt- 
kondensatoren verwendet, da sich in diesen am meisten « Joule pro Franken» 
aufspeichern lassen. Diese Kondensatoren kénnen an einem Schaltbrett wahl- 
weise in Serie oder parallel geschaltet werden. Da fiir unsere Messungen nur 
Impulsdauern von einigen Millisekunden und mehr in Frage kommen, sind 
wir nicht auf eine speziell induktionsarme Gestaltung der Zuleitungen zu den 
verschiedenen Kondensatoren und des Entladungskreises angewiesen (die 
Selbstinduktion unserer Hochfeldspulen ist von der Gréssenordnung 100 pHy 
bis 100 mHy). 


b) Das Aufladegerdat 


Ein Aufladegerat, das in der Lage ist, Stromstarken von 100 mA bis 500mA 
zu liefern, wird im allgemeinen fiir unsere Zwecke vollstandig gentigen, da in 
diesem Falle die Aufladezeiten bereits kiirzer als die entsprechenden Abkth- 
lungszeiten der Hochfeldspule werden. Unser Aufladegerat wurde so entworfen, 
dass die Aufladespannung vorgewahlt werden kann, um eine beliebige Ein- 
stellung der Héhe des Feldmaximums zu erméglichen. Wie man aus der 
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Schaltung von Figur 7 ersieht, wird hier der Ziindzeitpunkt des Thyratrons 
mit Hilfe der Steuerspannung am Gitter beeinflusst. Die Aufladespannung 
kann so zwischen Null und der Amplitude der sekundaren Wechselspannung 
des Transformators kontinuierlich eingestellt werden. 


Figur 7 


Prinzipielles Schaltschema der Aufladeapparatur. 


Wir sind in der Lage, mit diesem Aufladegerat unsere Kondensatorenbat- 
terie bis auf 2800 V aufzuladen, wobei Stromstarken bis zu 500 mA daraus 
entnommen werden kénnen. 


c) Der Entladeschalter 


Fiir unsere Zwecke kommen nur Schalter in Frage, die einerseits in der 
Lage sind, Stréme von der Gréssenordnung 1000 bis 10000 A zu schalten, und 
bei denen andererseits der Einschaltvorgang schon nach wenigen Mikrosekun- 
den abgeschlossen ist. 

Wenn vor allem auf die zweite Bedingung Gewicht gelegt werden muss, so 
wird sich im allgemeinen die Verwendung eines elektronischen Schalters 
(Thyratron oder Ignitron) empfehlen. 

Thyratrons und Ignitrons haben jedoch den Nachteil, dass sie sowohl die 
maximale Stromstirke als auch die Impulsdauer beschranken. Doch eignen sie 
sich besonders gut im Zusammenhang mit einer Elektrolytkondensatorenbatte- 
rie, da durch die Gleichrichterwirkung dieser Schalter eine schadliche, negative 
Aufladung der Elektrolytkondensatoren verhindert oder zumindest stark 
reduziert werden kann. 

Eine weitere Méglichkeit, das Schalterproblem zu lésen, besteht in der Ver- 
wendung von Funkenschaltern. Solche Schalter kénnen leicht selber hergestellt 
werden. Eine ausfiihrliche Beschreibung zweier solcher Schalter findet man in 
einer Arbeit von CULLINGTON, CHACE and MORGAN [11]. 
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Fiir unsere Apparatur haben wir ein Ignitron (P1 5555) alsEntladeschalter 
benutzt und damit sehr gute Erfahrungen gemacht. Dieses Ignitron soll nach 
den Angaben der Lieferfirma wahrend 0,15 s KurzschluBstréme bis zu 10000 A 
aushalten. Es arbeitet in einem Spannungsbereich zwischen 200 und 2000 V. 


d) Die Hochfeldspule 


Die Richtlinien zur energetisch giinstigsten Dimensionierung einer Hoch- 
feldspule sind bereits am Schluss des 2. Kapitels angegeben worden. 

Die Daten von drei Spulentypen, die hier besonders haufig verwendet 
werden, sind in Tabelle 4 zusammengestellt. Die Spule Nr. 1 ist fiir die Erzeu- 
gung von moglichst hohen Magnetfeldern vorgesehen, wobei allerdings eine fiir 


Tabelle 4 
Charakteristische Spulentypen 


a Draht H/I 


SESE a Gel Ca ade : Mitte 20009 S000 pF cine 
Win- , 2000 u : 2000 UF 
mm mm mm | dungen mm Oe/A Oe mHy ms 
iF 2.55 9,0 | 20,6 231 130 122 400000 0,18 On 
x 0525 
pe 6,6 254 fO15,5 7500 0,4 155 100000 420 30 
3 2 8,5 | 76,8 1085 1,30 176 230000 10,8 2 
x 025 


Die Spule Nr. 1 wurde mit zwei Glasfaserschichten verstarkt (siehe Seite 27). 


unsere Verhaltnisse sehr kurze Impulsdauer und ein nur in einem sehr kleinen . 
Volumen einigermassen homogenes Magnetfeld in Kauf genommen werden. 
Mit der Spule Nr. 2 wurde das Hochfeld im Gegensatz dazu wahrend einer 
relativ langen Impulsdauer und in einem grossen Volumen erzeugt. Diese Spule 
wurde so dimensioniert, dass sie tiber ein spezielles Glasdewargefass von 10 mm 
Aussendurchmesser geschoben werden kann. In diesem Dewargefass kann eine 
Messung zum Beispiel in fliissigem Helium (4,2° K) oder bei anderen Tempera- 
turen vorgenommen werden, wahrend sich die Hochfeldspule in fliissiger Luft 
befindet. Die Spule Nr. 3 endlich ist fiir Messungen vorgesehen, bei denen speziell 
Gewicht auf die Homogenitat des Magnetfeldes langs der Achse gelegt werden 
muss (Widerstandsanderung im longitudinalen Magnettfeld). 


e) Die Messung des Magnetfeldes 


Es hat sich gezeigt, dass fiir unsere Messungen das Magnetfeld mit genii- 
gender Genauigkeit der Stromstarke in der Hochfeldspule proportional gesetzt 
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werden darf, da Abweichungen von dieser Annahme, die aus den Deformationen 
der Hochfeldspule wahrend der Entladung resultieren, gréssenordnungsmassig — 
kleiner sind als die Fehler, die sich zum Beispiel aus der Nichtlinearitat und der 
Strahldicke des Kathodenstrahloszillographen ergeben. Um dies nachzuweisen, 
wurde fiir die verschiedenen Spulentypen jeweils die integrierte Spannung 
einer kleinen Pick-up-Spule gegen die Spannung an einem kleinen, induktions- 
freien Messwiderstand im Entladungskreis (Stromstarke), aufgetragen. Die 
Steilheit der Geraden, die so entsteht, liefert zugleich einen zuverlassigen Wert ~ 
fir den Quotienten H/I, denn die Berechnung dieses Quotienten aus der 
Geometrie unserer kleinen Hochfeldspulen wird im allgemeinen zu ungenau. 


6. Die Messung der magnetischen Widerstandsadnderung von Metallen — 


Unsere Messmethode arbeitet nach folgendem Prinzip: Durch die Metall- 
probe, die wegen der Erschiitterung der Hochfeldspule méglichst solid in dieser 
montiert werden muss, wird ein Gleichstrom geleitet. Der Spannungsabfall, der — 
dabei an der Probe entsteht, wird auf den Vertikalverstarker eines Kathoden- — 
strahloszillographen gegeben, wahrend zugleich die Spannung am Messwider- | 
stand des Entladungskreises an den Horizontalverstarker gelegt wird. Wahrend ~ 
einer Entladung wird so der Verlauf der magnetischen Widerstandsdnderung ~ 
direkt auf den Schirm des Kathodenstrahloszillographen gezeichnet (Figur 8). 


Bw 
SEererae 
Seer. 


: aH 
Hanieetieteste: 


a Figur 8 b 
Die Messung der magnetischen Widerstandsanderung von Cu (Q4 2/300 = 0,88%): 


@ Bei 77°K (Ay aq, = 260000 Ajom = 325000 Oe, A0/09 mag = 9,51). Die hochfrequenten Sté- : 
rungen werden durch die Erschiitterung der Hochfeldspule verursacht. 


b Bei 4,2° K (A, 44 = 160000 A/cm, A0/Oomaz = 5,82). Die Hysteresis ist auf Wirbelstromauf- 


heizung zuriickzufiihren. 


; Die durch das pulsierende Hochfeld im Messkreis induzierte Spannung 
dussert sich durch eine Hysteresis und wird hier mit einer Pick-up-Spule iiber 
ein Potentiometer auskompensiert. 

Mit dieser Methode sind hier zahlreiche Messungen der magnetischen Wider- 
standsaénderung von polykristallinen Metallen und Legierungen ausgefiihrt 
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worden [12}. Speziell wurde die longitudinale, magnetische Widerstandsaénderung 
von Ag, Au, Li, Ni, Pb, Pt, Sn wnd-Zn-bei 42°K untersucht [13], 
CHAMBERS [14] weist darauf hin, dass die ersten Messungen der magneti- 
schen Widerstandsinderung von Elektrolytkupfer, die hier in pulsierenden 
Hochfeldern bei 4,2°K ausgefiihrt wurden [15], nicht, wie es nach der Kohler- 
schen Regel zu erwarten ware, mit seinen (und Boroviks [16}), in Dauermagnet- 
feldern an hochreinen Cu-Proben vorgenommenen Messungen iibereinstimmen. 


“nsere neuen 
Messungen 
p= 88-10? 


zg 


FE 
BorOVIK R=217-10 8 3 7” “OLSEN und RINDERER 
T eee. ae f.-80-10 3 


: SS SS eee 
0 5 10 15 20 [Moe] 


Figur 9 
Transversale magnetische Widerstandsanderung Ag/g von Cu bei 4,2°K als Funktion von H/r7, 
wobei 77 = 04,9/0399- Die eingezeichneten Kreuze zeigen die Messfehler unserer Messungen an. 


Unsere Messungen (und auch eine Messung der magnetischen Widerstands- 
anderung von Al {17]) mussten aus energetischen Griinden mit sehr kleinen 
Hochfeldspulen und entsprechend kleinen Proben ausgefiihrt werden. Es 
zeigte sich denn auch, dass eine neue Messung an Elektrolytkupfer Resultate 
lieferte, die im Kohler-Diagramm gut mit den Messungen von CHAMBERS iiber- 
sinstimmen?). Die verschiedenen Resultate sind in Figur 9 zusammengestellt. 
Es scheint, dass nach den verschiedenen Verbesserungen, die im Laufe der 
letzten Jahre an unserer Messmethode vorgenommen wurden (gréssere Hoch- 
feldspulen, langere Impulszeiten, verbesserte Kompensationsmethode), nun 
zuverlissige Resultate erwartet werden diirfen. Unsere Messung wurde an 
siner Cu-Probe von 0,88°%, Restwiderstandsverhiltnis (@4,9/@399) bei 4,2°K in 
Magnetfeldern bis zu 160000 A/cm (200000 Oe) vorgenommen. 


4) Diese Resultate stimmen auch gut mit Messungen von DE Launay und WeBBeER tberein, 
n die uns freundlicherweise vor der Veréffentlichung Einblick gewahrt wurde. 
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Herrn Prof. Dr. P. GRASSMANN, Vorstand des Institutes fiir kalorische 
Apparate und Kaltetechnik, der mir auch mit einigen sehr niitzlichen Rat-— 
schlagen half, méchte ich an dieser Stelle meinen herzlichen Dank ausdriicken. 
Ausserdem bin ich den Herren Dr. J. L. OrsEN und B. Lirut fiir ihre Mitarbeit 
beim Entwickeln dieser Apparaturen und fiir ihre Hilfe beim Redigieren dieser 
Veroffentlichung zu Dank verpflichtet. Diese Arbeit wurde aus dem «Alumi- 
niumfonds Neuhausen» und aus dem «Arbeitsbeschaffungsfonds des Bundes» 
finanziert. Beiden Institutionen sei auch an dieser Stelle fiir die zur Verftigung ~ 
gestellten Mittel gedankt. 
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Summary 


The design of an apparatus for producing very high pulsed magnetic fields using 
condenser discharges is discussed. Efficiency, heating and strength are investigated. 
The calculations are compared with experimental data on coil types producing 
320000 A/em (400000 Oe) for 1 ms in a volume of 0,5cm*® and 80000 A/cm 
(100.000 Oe) for 100 ms in 10 cm’. Results of magnetoresistance measurements on 
copper to 260000 A/cm (325000 Oe) at liquid nitrogen temperature and to : 
160000 A/em (200000 Oe) at liquid helium temperature are given. In a Kohler 


diagram these results are in agreement with work by CHAMBERS on purer specimens 
in lower fields. 


(Eingegangen: 27. Juni 1959.) 
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The Plastic Indentation of a Semi-Infinite Solid 
by a Perfectly Rough Circular Punch 


By GrorGE Eason?) and Ricuarp T. SHIELp2), Newcastle upon Tyne and 
Providence, R. I. 


1. Introduction 


Although many problems of plastic flow under conditions of plane strain 
have been solved, there are, as yet, few solutions to axially symmetric problems. 
Previous work includes solutions to the problems of the partially plastic thick- 
walled tube (see for example [6]*), and the problem of converging plastic flow 
in a conical channel [9,7]. These problems have been solved for both the Tresca 
and the von Mises yield conditions. 

SHIELD [8] has considered the general problem of the plastic flow of metals 
under conditions of axial symmetry for the Tresca yield criterion and associated 
plastic stress-strain relations. The following work is essentially a continuation 
of that paper. In [8], as here, the material was assumed to be isotropic, non- 
hardening, rigid-plastic and a discussion was given of the types of plastic flow 
which could occur. The discussion and the applications in [8] showed that the 
plastic regimes which satisfy the hypothesis proposed by Haar and von 
KARMAN [4] are of importance in the solution of axially symmetric problems. 
This hypothesis states that the circumferential stress is equal to one of the 
principal stresses in the meridional planes during plastic deformation. The 
problems considered under this hypothesis were the incipient necking of a . 
circular cylinder stressed to yielding in tension, and the indentation of a semi- 
infinite solid by a smooth, circular flat-ended rigid punch. 

In this paper, section 2 contains a résumé of the equations of [8] relevant 
to the present work. In section 3 a possible velocity field is obtained for the 
compression of a circular cylinder. In sections 4 and 5 the indentation of a 
semi-infinite solid by a perfectly rough, circular flat-ended rigid punch is 
considered. A complete solution to this problem is obtained by numerical 
methods, the numerical work being performed on the high speed digital 
computers AMOS and FERDINAND. The value of the average pressure over 
the surface of the punch is found to be 6-05 k, where k is the maximum shearing 
stress permissible in the body, in contrast to the value of 5-69 & obtained in [8] 


1) King’s College, Newcastle upon Tyne. 
2) Brown University, Providence, R.I., USA. 
3) Numbers in brackets refer to References, page 42. 
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for the smooth punch. The maximum numerical value of the ratio between the 
shearing stress t,, on the surface of the punch to the normal stress o, is found 
to be 0-139. It follows that the condition of perfect roughness between the 
punch and the material can be replaced by the condition that the coefficient of 
friction between the punch and the material must exceed 0-139 for the solution — 


to apply. 
2. The Basic Equations 


The basic equations of axially symmetric plastic flow were discussed in — 
detail in [8] for a rigid-plastic nonhardening material obeying TRESCA’s maxi- 
mum shear stress yield criterion. As in the application of the theory in [8], the 
hypothesis of HAAR and von KARMAN [4] is assumed for plastic states, 1. e. the © 
circumferential stress is taken to be equal to one of the other two principal - 
stresses in the meridional planes. More specifically, the circumferential stress 
is taken equal to the algebraically greater of the other two principal stresses, 
as in the regions of deformation to be considered the radial velocity is positive 
and produces a tensile circumferential strain rate (plastic regime F of [8]). The 
Tresca yield criterion then requires the principal stresses in the meridional 
planes to differ by 2 k, where k is the maximum shearing stress permissible in . 
the material. A short summary is now given of the relevant equations. : 

In a cylindrical polar co-ordinate system (7, 0, z), the only non-vanishing ~ 
stress components in a stress distribution which is axially symmetric about the - 
Z-axis are 6,, 6g, G., T,,. The circumferential stress, o, is a principal stress, the 
other principal stresses 61, 6, (0, = 05) in the meridional planes being given by 


(1) 


For the plastic regime under discussion we have 
09 = 0,= 0,428. (2) 


The equations of equilibrium take the forms 


in the absence of body forces. 

Equations (2), (3) provide four equations for the determination of the four 
stress components so that in a sense the problem is ‘statically determinate’. 
The system of equations is hyperbolic with orthogonal characteristics in the 


Vol. XI, 1960 The Plastic Indentation of a Semi-Infinite Solid 35 


(r, z)-plane which coincide with the slip-lines, the lines of maximum shearing 
stress. The slip-lines will subsequently-be referred to as «- and f-lines. 

Surface elements perpendicular to the «- and f-lines are acted upon by a 
shearing stress & and a normal stress which will be denoted by —#. The stress 
components are given by 


o.=—p—ksin2p, o,=—pt+ksin2y, t,,=kcos2y, o,=—p+kh,. (4) 


where ¢ is the inclination of the «-line to the r-axis. The equilibrium equations 
(3) determine # and q, and referred to the slip-lines the equations take the 
forms r 
ds 
” 


dp + 2k dy + k (sing + cosy) —* = 0 on an eline , 


(5) 
dp — 2k dg — k (sing + cosq) sis =0on a f-line, 
where ds, ds, are elements of length along the «- and f-lines. 
On the axis of symmetry the shearing stress is zero and the radial and 


circumferential stresses are of equal magnitude, so that on the axis 
G=7m. (6) 


Using this result, an expansion of # and @ for small r shows that in the neigh- 
bourhood of the axis of symmetry 


p+4kw=const onan eline, p—4kgy=const ona feline , (7) 


to first order. , 

The velocity field is assumed to have radial symmetry, although this does 
not necessarily follow from the symmetry of the stress field. If u, w are the » 
velocities in the 7, z directions respectively, the non-zero strain rates are given by 


Ou uU Ow _ Ou Ow 
Sot igure Ample wele Geter TG ey (8) 
with ¢, a principal strain rate. The other principal strain rates €), & associated 
with the oj, o, directions are 
1 1 an 
Cia <Q (é, a 8.) +9 2Q { (é, ee EJP he et te 
; , 3 (9) 
ope 2 (é, oF é,] a { (é, ae é.) is Vis ia 
The incompressibility of the material requires, with (8) for the strain rates, 


Ou uU a 10 
RASA Oz 0. 0) 
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This equation and the isotropy condition, which can be written 


Ou ow 

pti re BER 11 
ga cot2@, (11) | 
rons 


serve to determine the velocity components #,w when the slip-line field is — 
known. For numerical work it is convenient to introduce the components of © 
velocity U, W along the «- and f-lines respectively, given by | 


5 


U=ucosp+wsing, W=—u sing + w cos@. (12) % 
Equations (10) and (11) on w, w require U, W to satisfy the equations 


dU — W dp + “S* =0 on anocline , | 
(13) | 
u dsp 


dw + Udp + 5- =0 ona fline. 


: 
For the velocity field to be associated with the stress field the velocities 
must be such that | 

: 


Ts See eee (14) 
i& 
where : 
_.0U , ow dy dp 
LS ggg pga ie des ears ey 


The theory outlined in this section will now be applied to obtain a defor- 
mation mode for the compression of a circular cylinder and to obtain the solu- 
tion to the problem of indentation by a rough circular punch. 


OD Compression of a Circular Cylinder 


In section 3 of [8] several examples were given of possible incipient defor- 
mation modes for a circular cylinder of rigid-plastic material stressed to yielding 
in uniaxial compression. The stress component o, then has the value — 2 k, all 
other stress components being zero. In this case the characteristics are straight : 
lines inclined at an angle of 2/4 to the co-ordinate axes. The cylinder occupies 
the region 7 <1, z = 0, and velocity solutions were obtained which are such. 
that the normal velocity over the end of the cylinder, z = 0, is constant. Here : 
we obtain a velocity field which satisfies the further restriction that the radial 
velocity w be zero over the end z = 0 of the cylinder. 

For simplicity the normal velocity is taken to be unity over the end OA of | 
the cylinder (Figure 1). The velocity field involves a rigid body motion of region 
OAC in the axial direction and flow in the region A BC. As in [8] it can be shown. 
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that w= w= 0 on BC because of the rigidity of the material above BC. The 
incompressibility condition and the-condition of zero shearing strain rate 
together with the continuity condition across AC then require «and w to be 
given by 


l 1 
get (i= tan*y)*, w= a cos“!(tany) , (16) 
in region ABC, 
u 1 
GOL Ce y = me Ws 


where the value of the inverse cosine lies between 0 and z. The angle yp is 
defined as in Figure 1. The velocity field defined by (16) was given previously 
in a different form in [8], apart from a multiplicative constant (equations 
(3.3) of [8}]). 


W=1 
Figure 1 
Cylinder stressed to yielding in compression, 


4. Circular Punch: Stress Field 


In [8] a complete solution to the problem of the indentation of a semi- 
infinite solid of rigid-plastic material by a smooth, circular flat-ended rigid 
punch was found. The closely allied problem of the indentation of a semi- 
infinite solid of rigid plastic material by a perfectly rough, circular flat-ended 
rigid punch is investigated here. These problems are both of importance in 
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testing the hardness: of metals (see, for example [10)). The body is assumed 
to occupy the region z = 0, with the origin of co-ordinates at the centre O of 
the punch (Figure 2). The punch indents the portion OA of the surface of the 
solid, z = 0. The plastic stress field near the punch is assumed to satisly (2) 
This assumption is shown to be correct by the fact that the solution to the 
problem is obtained. 


HTK Se SD, 
SQLS 
ZS 


0 | 
a 9 oto 
Som 2°90 Ss 8 ix} 
Ss ww Oo w+ 6 
S3xsz5 8 &§ $ SHS 
26 G6H G6 G6 GB HHH 
= 
5S 

Figure 2 


Slip-line field due to indentation by a rough circular punch. 


Referring to Figure 2, the stress field in region ABC is determined by the 
stress-free surface AB. The solution for this type of region has been tabulated 
in [8], and the results of Table I of that paper may be applied. The line AC is 
an «-line and BC is a f-line. The field in region ACD is determined by the 
g-line AC and the singular point A, the fan being terminated by the «-line AD 
which meets the axis of symmetry at an angle of 7/4, at D. The field in region _ 
ADO is determined by the «-line AD and the condition that the f-lines of the - 
region meet the axis of symmetry at angles of 7/4. The slip lines must meet the 
axis of symmetry at angles of z/4 in order that the shearing stress be zero on - 
the axis. 

The stress field in region ACDO was obtained by numerical integration of 
(5). These differential equations were replaced by finite difference equations 
and the same numerical procedure used as in [8]. Equations (7) were used to 
determine the position of the point D. Figure 2 shows the main details of the 
characteristic network obtained as a result of the numerical integration. The 
calculated pressure distribution over the punch is shown graphically in Figure3. 
The average pressure over the punch is found to be 6-05 k. For comparison, the 
pressure distribution found in [8] for a smooth punch is also shown in Figure 3, 
the average pressure over the punch in this case is found to be 5:69 k. With the 
distance OA as the unit of length, the distance OB, in Figure 2, was found to 
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be 1-88, and the distance OD was found to be 0:57, The angle of the fan was 
found to be 116°. A coefficient of friction may be defined along OA by writing 


Ty z 


poscre tas, (17) 
The maximum value of « along OA was found to occur at A and to have the 
value 0-139. Thus, the solution obtained here holds for a rough punch provided 
that the coefficient of friction between the material and the punch exceeds 
0-139. 


8-0 


7 


Figure 3 


Pressure distribution over punch. 


As in the case of the smooth punch the stress field has been extended into © 
the rest of the body. The method of extension closely follows that used in [8} 
in which theorems developed by BrsHop [1] for the case of plane strain were 
applied. 

The extension of the stress field is indicated in Figure 4, and was obtained 
by assuming the material to be fully plastic and satisfying (2). The f-line BCD, 
together with the condition that slip-lines meet the axis of symmetry at angles 
of 2/4, determines the field to the left of the «-line through B. The field to the 
right of the a-line through B is determined by the «-line through B and termi- 
nating the field by a stress-free surface, BJ K. The a-lines begin to intersect one 
another at G, a point on the «-line through B, and a shock is introduced, shown 
by the heavy line in Figure 4. Another shock is introduced at point M, where the 
x-lines intersect on the «-line through C. The two shocks intersect at point H, 
and the stress field was continued by the introduction of a shock HL, and a 
fan HJ of small angle (shown by the broken line in Figure 4). The stress-free 
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boundary is parallel to the axis of symmetry at the point K. The stress field 
was terminated by the introduction of a stress discontinuity KLN , above which 
the stress is uniaxial compression or tension parallel to the axis of symmetry. 
The stresses in the region are indicated in the figure, and do not violate yield. 
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Figure 4 


extension of the slip-line field into the rigid region. 


The total force acting on KLN was calculated and found to be within 3-5°% 
of the pressure on OA. This discrepancy is almost certainly due to the inaccu- 
racy of the finite difference forms of (5) near the axis, and of (7) away from the 
axis. Neither set of equations is really satisfactory in the neighbourhood of the 
axis of symmetry. With the distance OA as the unit of length, the point K is 
distant 3-59 from ON, 2:21 from OB, and N is distant 3:65 from OB. 

The stress field obtained is statically admissible and the value of 6-05 2 is 


a lower bound for the average pressure by the limit analysis theorems due to 
DRUCKER, GREENBERG and PRAGER [3]. : 


5. Circular Punch: Velocity Field 


An incipient plastic velocity field is obtained by assuming that flow id 
confined to the region A BCD of Figure 2. It is assumed that the punch has unit 
velocity parallel to the axis of symmetry, and that the region OAD moves as a 
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rigid body in the axial direction with the velocity of the punch. The velocity 
field in region ABCD must satisfy the conditions that the normal velocity 
across AD is continuous and that the normal velocity across BCD is zero, as 
the region above this line is assumed to be rigid. As AD is an «-line, the first 
of these conditions implies W continuous across AD, so that 


W = cosp on AD. (18) 
Along AD the first of equations (13) holds, and with (18) this takes the form 
dU — cosp dp + (U —sing)  =0. (19) 

The general solution of this equation is 
urane (20) 


where A is a constant. To avoid an infinity in the value of U at the point D, 
A must be taken to be zero, so that 


U =sing on AD. (21) 


Thus the velocity field is seen to be continuous across 4D. A similar argument 
to the above shows that 


U=W=0 along BCD, (22) 


so that the velocity is also continuous across BCD. 

The point D is a singular point of the velocity field. Near this point the 
slip-lines are inclined at angles of 2/4 to the co-ordinate axes, and the velocity 
field given by (16) may be taken to hold. This velocity field was assumed to 
apply in the region DEF of Figure 2. Since the velocity is now known along 
the lines AFECB, the velocities within this region may be determined by 
integration of (13). The point A is also a singular point of the velocity field. 

The velocity field in region AFECB was obtained by numerical integration 
of (13). The same numerical procedure was adopted as in [8], the differential 
equations being replaced by finite difference equations. The incompressibility 
of the material provides a check on the accuracy of the calculation. It was 
found that the flow across AB was within 1% of the flow across OA. The 
deformation of an initially square grid which would result if this incipient 
velocity field were maintained for a short interval of time is shown in Figure 5. 

The velocity field must also satisfy the inequalities (14) if it is to be asso- 
~ ciated with the plastic stress field. The inequalities were checked numerically 
at various points of the field, a sufficient number of points being taken to 
ensure that the inequalities held everywhere. 
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The velocity field is a kinematically admissible deformation mode and 
application of the limit analysis theorems shows that the value 6-05 # is an 
upper bound for the indentation pressure. The value 6-05 & is thus the actual 
indentation pressure, as it is also a lower bound. A theorem due to Hirt [5] 
states that where deformation is actually occurring the stress field is unique. 
Also a theorem to BisHop, GREEN and Hirt [2] states that if any region of 
a complete solution is necessarily rigid, then it must be rigid in all complete 
solutions. Application of these theorems shows that the stress field of Figure 2 
in region A BCD is the actual plastic stress field, and plastic deformation can 
only occur in region ABCD of Figure 5. 


Axis of symmetry 


i} 
| 
ae See 


Figure 5 


Resulting deformation of a square grid if the incipient velocity field were maintained for a short 
period of time. 
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Zusammenfassung 


In dieser Arbeit werden das plastische Spannungsfeld und ein zulassiges Geschwin- 
digkeitsfeld fiir den eben begrenzten Halbraum gegeben, der unter dem Einfluss 
eines ideal rauhen, starren Stempels mit kreisférmigem Querschnitt steht. Das 
Material ist als starr-plastisch vorausgesetzt, ohne Verfestigung, und der Fliess- 
bedingung von TreEsca geniigend. Es wird gezeigt, dass die Hypothese von Haar 
und von KARMAN auf dieses Problem anwendbar ist, wonach zwei von den drei 
Hauptspannungen gleich sind. Es wird auch eine giiltige Fortsetzung des plasti- 
schen Spannungsfeldes ins starre Gebiet in der Nahe des Stempels erhalten. 
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Flow of an Electrically Conducting Fluid Past a Porous Flat 
Plate in the Presence of a Transverse Magnetic Field 


By ANADI SHANKAR Gupta, Kharagpur, India!) 


1. Introduction 


In recent years, the subject of hydromagnetics has attracted the attention 
of many authors in view not only of its own interest but also of its enormous 
applications to problems of geophysical and astrophysical significance. 

The general equations of unified velocity and magnetic fields for the flow 
of an incompressible, viscous and electrically conducting fluid subject to a 
magnetic field have been derived by BATCHELOR [1]?). 

In this paper, we have discussed the effect of a transverse magnetic field 
on the steady flow of an incompressible electrically conducting fluid past an 
infinite porous flat plate when the plate is subjected to either suction or injec- 
tion. Exact solution has been obtained for the modified Navier-Stokes and 
Maxwell equations under the usual assumptions of magneto-hydrodynamics. 
The equation of heat transfer including viscous and Joule dissipation has also 
been integrated and the rate of heat transfer calculated. The corresponding 
problem in the absence of magnetic field has been solved by GrirritH and 
MEREDITH [2]. 


2. Basic Equations 


The fundamental equations of hydromagnetics are: 


(a) MAXWELL’s equations: 


1) Indian Institute of Technology. 
2) Numbers in brackets refer to References, page 90. 
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divB=0, (2) | 
0B 
curl B= = ——, “hyd 
&.  0E 
eid peas ete 4 
curl = 40 je ene) (4) 
B= un, H. (5) q 
(b) Oum’s law: 
j=o(E+q x B)+0q. (6) 


(c) The equation of continuity for an incompressible fluid : 
a iv == Oe (7) 
(d) The modified Navier-Stokes equation: 


$4 4 (q+ grad) q= — ; grad +9 V2q + a] x B+ £ E. 8) 
Here, Eis the electric field, B the electromagnetic induction, H the magnetic 
field, j the current density, Q the excess charge density, e the dielectric constant, 
u, the magnetic permeability, o the electrical conductivity and c the velocity 
of light. All electromagnetic quantities are measured in electromagnetic units. 
The quantities e and yw, are assumed to be uniform throughout the field. 
ELSASSER [3] has pointed out that the ratio of the displacement current to 
the conduction current in a large class of phenomena is so small that displace- 
ment currents can be neglected entirely. In our case it is zero simply by 
virtue of the fact that the motion is steady. Consistent with the above assump- 
tion, the convection current Q q in (6) and the electric force Q E/o due to excess 
charge in (8) are also ignored. In the following analysis the equations (4), (5) 
and (8) should be understood to have been simplified by the above assumptions. 
Eliminating E between (3) and (6) with the help of (2), (4) and (5), we have 


OH 
@ Ul x A)+nV A, (9) 


where 


ee (10) 


40,0’ 


which is interpreted as a magnetic diffusivity. Again elimination of j between 
(4) and (8) gives 


0q 1 
ae + (4 grad) q = PKs. grad p + » V7? g + ae (curl x) een 
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Equations (9) and (11) together with (2) (i. e. div H = 0) and (7) constitute the 
basic equations of hydromagnetics. 

In the present problem, we take the coordinate origin O at any point on the 
infinite plate, which is assumed to be nonmagnetic and electrically non- 
conducting and use Cartesian co-ordinates x and y with the axis OX along 
the plate OY and perpendicular to it. The liquid flows past the porous plate 
along OX and a transverse magnetic field H, is applied to it. 

For steady flow past an infinite plate, we can assume that all physical 
quantities are functions of y only. Let « and v denote the components of the 
vector q in the directions of the co-ordinates X and Y. 

The equation of continuity then becomes 


dv 
wer 0:, (12) 
which gives 
) =v, = const. (13) 


The modified Navier-Stokes equation (11) then gives 
du au a dH 


0° Gy Hat * an OY dy? ie 
_ ap lle dH, i 

4 dy 4% i. dy. ~ 2) 

Again the magnetic equation (9) gives 

: aH 

jy UH, — %H,) + ae =0, (16) 

aH 
Ragan (17) 

The boundary conditions for the velocity field are 

y= OF 4 = 0, v= 5; yoo: u=U, (18) 


where U is the velocity at far distance from the plate. 
The boundary conditions for the magnetic field are that both the tangential 
component of the magnetic field and the normal component of the magnetic 
induction are continuous at the plate. Since, however, both the fluid and the 
: plate are assumed to be nonmagnetic in the present problem, the continuity 
_ condition for the magnetic induction may be replaced by that for the magnetic 
field itself. The total magnetic field at the plate is the vector sum of the 
imposed and induced fields, the latter resulting from the induced current 
distribution in the fluid. This distribution is obtained from the solution and is 
therefore unknown. On the other hand, since the velocity tends to the constant 


46 ANADI SHANKAR GUPTA ZAMP — 


value U asymptotically at far distance from the plate, no interactions exist | 
there. Hence the magnetic field tends to the imposed field H, as y > oo. The 
boundary conditions for the magnetic field may therefore be taken as 


y=0: H,=Hy; yoo: Bye gwigs Us (19) 
These imply that the total flux of current in the direction OY is zero since 


4 il 
fae curlH. 


3. Integration of the Equations 

The solution of (17) subject to (19) is 
H, = const = Hy, (20) 3 
while (16) with H, = Hy yields on integration 


“4H,—vH,+n a = const = U"HS ; (21) | 


where the constant is determined from the fact that H, >0, dH,/dy >0 
and u—> U asy> oo. 
Again (14) with (20) gives on integration 


du 
BRA t he 5 +42 HoH, + o%U; 


where the constant of integration is evaluated from the fact that u >U, 
du/dy > 0 and H, >0 as y + oo. The above equation gives 


H,= 


z ae E My (w— U)— ps al (22) 


which gives the induced field when the velocity field is known. Substitution 
of H, from (22) in (21) gives 


au Uo du lt, He ue H, U (4 & \ 
aE a oe late Ae | Me Ho 0 v2 | 
dy? mp8 H) dy lana 7 | 4a \ m, Hy Ho). a 


If the pressure at infinity is p,,, it follows from (15) that 


Me 2 
p a 87 He i Ps , 
which shows the balance between the fluid and ‘magnetic pressure’ because 
H, has been shown to be uniform throughout the whole field. 
We shall now consider two cases. 


Vol. XI, 1960 Flow of an Electrically Conducting Fluid’ Past a Porous Flat Plate 47 


Case F? v4, <0 


This case corresponds to suction at the plate. Putting vy = — v (v > 0), the 
solution of (23) is given by 
n= U+ Aexp{s” Ey 1(no+p) + or (no+p ote ea | 


+ Bexp eg eed ome (yo +p) — ex (no+p)?+4 ee _ 0 0? | nu(| : 
Clearly the expression in (24) is real for 
2 (Me Hy 
wv (no + p+ 4 fae —o v*) 1 bs 


can be written as v? (70 — m“)?+ uw, Hj/x, which is positive. Again when 
Lt, Ho/4 x > o v? the expression 


Barer cena cs ee 4 
Yer ot mrt (MBE — 9 vt) ny 


becomes greater than v (7 9 + mu). Hence in this case the solution (24) subject 
to (18) gives 


u } (9 p\Pir , FR 02 
et wife Q 0 
parole hess Go ee]. es 


It is convenient to define some non-dimensional parameters such as 


local Reynolds number: Re = “2 : 
Ma f . — Me Ay : 
gnetic pressure number: S = Tae (26) 
together with the non-dimensional number Ry = »/n. 
Hence (25) can be written as 
u vy a 
Ces exp | ie [1 +R +VG—RPe+4S Rill, (27) 


which reduces to the ordinary asymptotic suction profile u/U = 1—e~*%?” of 
GrirFitH and MEREDITH when S = 0, i.e. Hy = 0. The velocity distribution 
given by (27) is plotted for different values of S when Ry = 0-5. It is seen from 
the graph that the effect of increase of S, i. e. the magnetic field is to increase 
the velocity. The point which seems to be of some physical interest is that 
the magnetic field should be strong enough to satisfy uw, 5/40 > ov? (ie. 
the magnetic pressure number S based on suction velocity should be > 1) in 


order that an asymptotic solution of the type (26) is possible. 
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If, however, u, H3/4a <ov%, i.e. S<1, then both the powers of e in 
(24) become negative so that the constants A and B can not be determined 
from the boundary conditions (18). Thus in this case the problem becomes 
indeterminate. The problem, however, can be made determinate by suitably 
changing the boundary conditions viz., by considering two parallel plates. 

The skinfriction t,, at the plate is given by 


re=M (Sug = ay + Rot VO-RPF4ASR], 8) 


which increases as S increases. 


Case liste Y 


This case corresponds to injection at the plate. Here again for 
u, Ho/4 a > o vg the solution of (23) subject to (18) is 


™ — 
wal exp [32 E (no+")- /8 ne — pet He null. 


Introducing S and R, this can be written as 


oy =l—exp{— 2% [V0 = Ry)? +45 Ry — (1 + Ro] (29) 


for S > 1. This shows that w increases with S. . 
The solution, however, remains indeterminate for S < 1 as in the previous 
case. 


Using (29), the induced field component H, can now be calculated from 
(22) as 


= a xP (— FW = Rt + 4S Ry (1+ R)}) | 
(30) 


x [1+ 7 WOR +45R,— (01+ Rl]. 


The value at the plate is given by 


j= =e 1 py pee = Res oR, See 


Using the definition of S the above equation can be written as 


iL 
SS eS awe / = 
EM) a tt yg VA- RP +4SR,— (1+ R)} 
410 


io 3 =f(Ro S), (31) 
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say. Since 


I 


2 j(Ry, §) = — = ROVE RP EASR, + (1 ~ Ro) 


for Ry < 1, it follows that the induced field decreases with increase in the 
magnetic pressure number. Since S = mu, Hj/4@v}, decrease in v, for a 
given H, will result in increase of Sy. Hence it follows that the induced field at 
the plate decreases with decrease of the injection velocity vp. 


4. Integration of the Equation of Heat Transfer 
The equation of heat transfer including viscous and Joule dissipation (for 
finite electrical conductivity o) is 
dT 


@T  j du\2 
ee sae A ay em, (= , 


where c, is the specific heat at constant volume and A the thermal conductivity. 
Since 


A 1 
y | aaa Fie curl H 


the above equation can be written as 


rig ® a2T al aH, 2 du\2 
0% GA or +aeats (et) te(Z)- G4) 


Substituting the values of w and H, from (29) and (30) in (32) we get 


Sas. tet at big MU? VE pe) ,-(2r6/»)v 
ay? Dae dy ey ( 16 w20 y2 # ‘ ; 3a) 
where 
1 = 4naovju fp (1+ B) 
p=+[VO— RP T4SR,— (1+ Ri], yo APE | 34) 


The solution of (33) subject to T = T, at y= 0 and T finite at y = oo is given by 


y? + 16 a? wp U? v3 Bo fee e~ (2%oBir) y) : (35) 


L—To= 32 m2 o v2 (2 B21 + cy m B) 


For o = «, the rate of heat transfer at the plate is given by 


Since for an incompressible fluid, c, =, it follows from the above equation 


ZAMP XI/4 
a 
‘g 


> 
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that 


Pe a (36) 
QU'a, — Zp+ Pr 


(where Py is the Prandtl number y ¢,/A), which increases as # increases. 
I should like to thank Prof. B. R. SETH and Dr. G. BANDYOPADHYAY for 


their valuable comments. 
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Velocity distribution over the plate for different values of S when Ry = 0°5. 
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Zusammenfassung 


Fiir die Strémung einer inkompressiblen, zahen, elektrisch leitenden Fliissigkeit 
liber eine pordse Platte werden die Gleichungen von NAVIER und STOKES, von Max- 
WELL und die Warmeiibertragungsgleichung integriert. Dabei wird vorausgesetzt, 
dass Fliissigkeit aus der pordsen Platte austritt oder durch sie abgesogen wird. Man 
findet, dass sowohl die Geschwindigkeiten als auch die Schubspannungen mit stei- 
gender magnetischer Druckzahl fiir Absaugen und Ausblasen anwachsen. Das 
induzierte Feld an der Plattenoberflache wird kleiner, wenn die Ausblasegeschwin- 
digkeit abnimmt. Fiir den Grenzfall unendlicher elektrischer Leitfahigkeit wachst 
der Warmeiibergang mit der magnetischen Druckzahl im Falle des Ausblasens. 


(Received: May 20, 1959). 
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The Possibility of Normal Shock Waves on a Body 
with Convex Surfaces in Inviscid Transonic Flow 


By GeorceE E. Gapp, Teddington, Mddx, Great Britain!) 


Introduction 


In inviscid transonic flow past a body with a convex surface, it can be shown 
by an argument essentially the same as that of Z1EREP [1]?), that shock waves 
normal to the surface cannot occur if the Mach number upstream of the shock is 
close to 1, and the flow is free from singularities. It would, however, be surprising 
if such shocks could not, in fact, exist in inviscid flow, in view of Emmons’ relaxation 
solutions [2, 3} which require them. Emmons suggested that the flow at the surface 
immediately downstream of the shock must have a singularity of the same kind as 
occurs for low-speed flow over a surface with a discontinuity in curvature. Such a 
flow has been considered by KoOpPENFELS [4]. The present note shows, ‘in confir- 
mation of Emmons’ suggestion, that a simple modification of KOPPENFELS’ solution 
appears to be applicable immediately downstream of a normal shock on a convex 
surface, making the occurence of such shocks possible. 


Analysis 


It can be shown (see for example [5]) that the deflection angle ¢« of the flow 
passing through a shock which is inclined at an angle 9 to the normal to the up- 
stream flow, as in Figure 1, is given by 


Y Pi 019% M3 q : 2 Pi aly 
oe ar 1 : Pr» = Pr s (Pe Fa ps)” 2 (Po > Pi) : 
Zy P1 Pia di 017i 


i ve 


(be =P [ary —1 — 24D Pe) 
(1) 


sin? ¢ = 


where suffices 1 and 2 denote conditions upstream and downstream of the shock, 
respectively, p is pressure, g is density, y is the ratio of the specific heat at constant 
pressure to that at constant volume, M is Mach number, and g is speed. Also 


ies cadet? (2) 
Po — Pi 


The denominator of the right hand side of (1) is approximately equal to 1 if the 


tan 0 = 


; 1) National Physical Laboratory, Department of Scientific and Industrial Research, Aero- 
_ dynamics Division. 
2) Numbers in brackets refer to References, page 58. 
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upstrearn Mach number is not much above 1. Hence the expression 


(y + 1) (b2 — Pi) (3) 


S = M2—1 
: 2 pi 


cannot become negative if sine is to be real. | 


Flow 
ee 


\ 
\ 
Suffix 1 conditions ' 


“ 


~_ Suffix 2 conditions 


Suffix a conditions 


Figure | Figure 2 
Flow through an oblique shock wave. Shock on a curved surface. 


Suppose the flow takes place over a curved boundary, and that the shock 
intersects the surface at a point O as in Figure 2. Let O be the origin of co-ordinates 
x, y parallel and perpendicular respectively to the surface at O. Denote by suffices 
a and b respectively conditions at the surface immediately upstream and down- 
stream of the shock. Then 


mel see 
(1 + HS Mi) - 1, (1+ 45 ~ M3), 


ee 


where J is enthalpy. Hence since the upstream flow is isentropic, so that 
aL aes 
I, Pa , 


Ssyar lle a)” 50a ae 


At O, S must be zero since the flow at the surface cannot abruptly change direction 


and hence sine must be 0. Since p, = ,, py = fp, at the point O, the condition — 
S = 0 requires 


equation (3) may be written 


bie ea (ME AY (5) 


This, of course, is the well-known relation for the pressure rise thr 
’ O h 
shock. From (4) it follows that at O E ugh a normal 


(Mi—1) dp 
iaceaa) ay pannel 


dy 2 Py 


ee 8) (=e _ ar» 
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or approximately 
GS\- (p+ Ty f/ dp, (dps \ \ 
(a ), ie Tye as) +a ly > 


if My is close to 1. This expression must not be negative if S is not to become 
negative for small values of y. 
Let x,(y) denote the x-coordinate of the shock wave. Thus if the shock is concave 


to the oncoming stream, as shown in Figure 2, ¥, is negative. Then 


B-() + (0) a 
and 
se — (+ (By ° 
From (7) 
(ad 0 
feats My ea °) 


because at 0, dxy,/dvy = 0, the shock being perpendicular to the boundary. However 
(dp,/dy), may not be equal to (0p/0y),, since (0p/dx), may be infinite and the product 
(0p/0x), (dx,/dy), may accordingly be non-zero. It is possible to relate (0p/0x), to 
the shape of the shock near O, as follows. 

The equations of momentum and continuity are 


Pe ala Ou Op 


ox 78” Oy =~ Ox? | i, 
v ae Cote, OF 
BO ii et By ay? (11) 
and 
0(o u) 0(@ v) 


where wu and wv are the velocity components in the x and y-directions. Since the 
Mach numbers upstream of the shock are not much above 1, f/g” is virtually 
constant downstream of the shock, as well as being exactly constant upstream of 
it. Hence the flow is approximately isentropic and irrotational everywhere, and 
for the downstream flow the continuity equation becomes at O 


Ge), - os (Ge), € al T os Ce ney 


Since UV» = Vit follows from the momentum equations that (0p/0x), = — OQ» Uy (Ou/Ox), 
and v, (0p/dv), = 0%), so that 

. (5) Pires (=) (13) 
i Ox }b (1 — Mg) \ ov Js 


od 
; 3) It can be verified that these relations are true both when +, y3/? for small values of y, as 
~ 

‘ZiEREP assumed would be the case for flow over a concave surface, and also when the downstream 


flow has a singularity of the type investigated by KopPENFELS and considered below. 
2 


ae “*N : 


54 Kurze Mitteilungen — Brief Reports — Communications bréves ZAMP 
and . 
OP:) ges COED (57) (14) 
Ox } b (1 — M@) \ oy Jo | 
Now 


av, 2/08 re (<2) dx, | 

payne (55 ), ox }2 dy | 

so that at O where dv,/dy = 0, | 
dy Jo Oy /b | 


Also, from Figures 1 and 2 it can be seen that v2 = uy (e + 6 — dx,/dy) and . 
0, = U, (6 — dx,/dy), so that vz = Uy V4/U, + M2 €, and 


[2 
dv,\ | Uy (=) 2 (Ma — 1) ee Ae 
(as ~ tg NOyla = +1) ME dy Jo (19) 


by equations (1), (3), and (5). Also since dx,/dy = 6 — 04/1, it follows from equations 
(1) to (3) that 


gE sy eis het ages (17) 


Hence from equations (6), (8), (9), (14), (15) and (16) it follows that approximately 


. 
: 
(45), = — ED 1/28) « (2)] — (sun dy 


since M2 — 1 ~ 1 — M3? and M} p, wy, ~ Mi p, u, if M, is not much above 1. Hence 


rae ersaibeatheste| 
ES eee é 18 
erg tee Seon thet Hp ee 
If the boundary past which the flow takes place is denoted by vy = h(x), then 


h(0) = h’(0) = 0 and in the vicinity of O, h(*) = h"(0) #2/2. Hence on the boundary 
close to O, v ~ h"(0) xu, ~ h"(0) x w, and from equation (11), 


— == 0, UE h’(0) ~ — 6, 02 bh" (0) 


It seems, therefore, at first sight that both (0p/0y), and (0p/dy), must be approxi- 
mately equal to — g, uj h"(0). Equation (18) thus gives approximately 


(=) 2 +10) 

dy }o was 

If the surface is convex to the stream, h”(0) is negative, so that (dS /dy), is negative, 
and this, as pointed out above, is incompatible with the shock relations. Thus it 
follows from (17) that for small values of y 


5) 


Bem + (FEE wr 9)™, (19) 


dy 


which is imaginary when the surface is convex. It might seem, therefore, that 
normal shocks cannot occur on convex curved surfaces. 
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However it is possible to envisage a singularity in the downstream flow of such 
a kind that (0p/0y), is not equal to — 0,-w>-h"(0) even though 0p/dv ~ — 0, u2 h" (0) 
when y = 0 and x is very small but non-zero. Such a flow would be similar ‘to the 
low-speed flow considered by Koppenretis [4], past a wall with a discontinuity 
im curvature, as in Figure 3. The velocity components close to the discontinuity 


Figure 3 


Boundary with a discontinuity in curvature. 


approximate, according to KoPpPENFELS’ solution, to the form 


=1—2(c, + ¢,) x log,c, (x2 + y*)2/24 2 [= Cy + (¢, + Cy) tan} | y (20) 
iy 5 
and 
= = 2 (c, + Cy) y log,cs (#2 + y?) 1/2 + 2 [ €y + (€, + cy) tan a ee (eye 
= 3 


where c,, ¢,, and c, are constants with the dimension (length)~?. It may easily be 
verified that these relations satisfy the condition for irrotational flow, 0u/dy = dv/dx, 
as well as the continuity condition du/dx = — dv/dy. For y = 0, v/uy = 226g 4, 
when x is positive, and v/u, = 2 2 (c, + 2c.) x when x is negative, whilst u ~ ™, 
for small positive or negative values of +. Hence (20) and (21) represent a flow 
past a wall with a radius of curvature equal to 9, = 1/2 a (¢, + 2¢,) for small 
negative values of x, and a radius 9, = 1/2 mc, for small positive x, as shown in 
Figure 3. 

This solution may readily be extended to a case where compressibility is im- 
portant. Consider the relations 


u 2 (cy + C9) ¥ ; 
ru, Ab ge: 5 = WEE log, cs [¥? oh (1 _ Mj) y*] = 


2) 1] (22) 
1. Me 2 
ae [= Cy + (C1 + 62) or a bes | y 
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= 2 (ey, + €) (1 — M3)4!? y log,c, (x? + (1 — M2) y7]1/? | : 
(23) 
a 
: 
These are of course related to the incompressible forms (20) and (21) by the ordinary | 
Prandtl-Glauert rule; this would not normally be valid for transonic flow with — 
M, close to 1, but it is so sufficiently close to the origin in the present case because 
of ae boundary conditions on v and because of the logarithmic terms. The aan 
tational cendition du/dy = dv/0x is still satisfied, whilst 0u/0% = —(dv/dy)/(1 — M). 
Thus condition (13) is satisfied, and it is therefore possible that relations (22) 
and (23) may represent conditions downstream of the shock. For y = 0, 
v/u, = 2c, x when »* is positive, as for the low-speed solution. Thus the radius 
of curvature of the surface downstream is equal to 9, = 1/2 a c¢,, where a positive 


value of g, means that the surface is concave to the stream as in Figure 3. If, then, 
the equation to the boundary is y = h(*), as before, ge, = 1/h"(0) od 


ewes y 


+2 [meq + (4, + cy) tan-t- 2 


t= 


: 
h"(0) 
Se (24) 
There are two other disposable constants, c, and c;, and of these c, may be found 
as follows. 
Conditions denoted by suffix } in the earlier analysis can now be considered to 


be the limiting values obtained by setting 7 = 0 in the appropriate quantities 
derived from equations (22) and (23), and then making y very small. Hence, since 
from equations (15) and (16), 


; Peer 1/ 
($7) ne 2 (M2 1) (4S -} 
0 


Oy }b (vy + 1) M2 dy 


it follows that 


dil (y +1) M2 
. dy ar, 5 (M2 nr i) (cy at C9) (1 a M?) 1/2 log, cs (1 ca. M3) 12 y 


for vanishingly small y. Thus S1/?2 « y log,cs (1 — M?2)1l2 4, and 
BS Wg. wit (2) 
dy 


as y > 0. Hence from equation (18), 


aah ei (a pak 


(22) = —ecutano 


But from equation (11) 


whilst 


(5°), = — 0 Ui (3 70,4 04) , 
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since downstream of the shock close to O 


me 1 — Me2)1l2. 
7 = 2 Oy Uy [= Cy + (cy + ¢,) tan ( M)) "| 


x 
by (11), (22), and (23). Hence, since 9, uw? ~ @, v2, 

3.%C,+ m6, = —h"(0) 
and from (24), 


cy = >> — h"(0). (25) 


Thus, if h"(0) is negative, as it is for a convex surface, the downstream flow is 
effectively the same as flow past a boundary of the shape shown in Figure 4, with 


Flow ly 


Radius Q,=-1/3h"(0) 


Figure 4 


Equivalent boundary for the downstream flow. 


@, = 1/2 a (¢, + 2¢,) = — 1/3h"(0), and — g, = — 1/h"(0). However the solution. 
only applies, of course, to the shock-wave case for positive values of «. 
The x-component pressure gradient downstream of O can easily be shown from 
(10), (22) and (23) to be 
Op 2 OM 


Ox a (1 aa M2) 1/2 a “ig Co) log, cs [4 an (1 a M75) pele : 


Hence along the surface 


Op 4 0, UB p : 
; eo, x 26 
Ox x(t — Maile © (0) 1086s 9) 


by (24) and (25). If the surface is convex (26) shows that the pressure diminishes 
very rapidly along the surface downstream of the shock, although at the shock 
the full normal-shock pressure rise is attained. This is in line with Emmons’ rela- 
-xation solutions. : 
Thus the argument that normal shocks cannot occur on a convex surtace 1s 
shown to be inconclusive. It has still not been established of course that they can 
exist, though Emmons’ work would suggest this. If they can exist, it seems quite 
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likely that relations (22) to (26) would be valid not only for convex surfaces, where 


h"(0) <0, but also for concave surfaces. In other words equation (19) may be ~ 


invalid even for concave surfaces. A more elaborate investigation, considering the 
second derivatives of quantities along the shock wave, would be necessary to 
settle this point. 
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Zusammenfassung 


Es zeigt sich, dass in der Tat senkrechte Stosswellen an einem Korper mit 
konvexen Oberflachen bei reibungsfreier transsonischer Stromung auftreten kénnen, 
wenn diese Strémung unmittelbar stromabwarts des Stosses eine Singularitat 
besitzt und diese zur Singularitat in der Stromung iiber eine Grenzflache mit Kriim- 
mungsdiskontinuitaét ahnlich ist. 
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Numerical Solutions of the Van der Pol Equation 


By Was-ey S. KroGpaut, Lexington, Kentucky, USA?) 


It is well known that the Van der Pol equation 
H+HKH=pF# (1 — x) 


has limit cycle solutions for all 4 => 0 [3, 4]?). The period, amplitude, and certain 
other parameters of the limit cycle solutions are known in asymptotic form in the 
limit of very small or very large yu [5, 6]. However, the characteristics of the solu- 
tions for intermediate values of « are to be had only by obtaining the limit cycle 
solutions numerically. The desirability of obtaining them has been pointed out by 
LEFSCHETzZ [2]. 

The solutions for integral « from 1 to 10 inclusive have been obtained by 
numerical integration using GILL’s modification [1] of the Runge-Kutta method. 
The computations were carried out on the IBM 650 calculator of the Northwestern 
University Computing Center. Subsequent independently programmed check cal- 
culations were made on the 650 calculator of the University of Kentucky Computing 
Center. 

The limit cycle in each case was identified as that solution whose half-period 
amplitude was equal to the initial amplitude. The half-period was determined as the 
value of the independent variable ¢ at the first zero. of # = dx/dt, the integration 
having been started in each instance at a trial amplitude +(0) = a), where #(0) = 0. 

The results are contained in the three graphs and Table 1. The first graph is a 
plot of x as a function of f¢, the different curves corresponding to the indicated values 
of w. Similarly, the second graph shows # against ¢, the ‘velocity’ as a function of 
‘time’. The third graph is the phase plane diagram of each solution over a half cycle. 

Table 1 compares the computed amplitudes a) and semi-periods P,/2 with the 
same quantities calculated from the asymptotic formulae for large u. To obtain 
some notion of the accuracy of the computed values, certain additional calculations 
were made. Thus, sample calculations were re-run using a double integration inter- 
val. Negligible differences in the two solutions indicate that the original interval 
was sufficiently small that truncation errors are of no consequence. On the other 
hand, by transforming the independent variable to t = w#, we are assured of negli- 
gible truncation errors but, because of the w-fold increase in number of steps, an 
estimate of the round-off error becomes possible. Only a coincidental compensation 
of the two sources of error could lead to the same solutions, and that for probably 
no more than one value of #. Such considerations offer reason to think the computed 
values correct to as many figures as given. 

The formulae used for the asymptotic values are 


16 


a*® = 2-0 + 0-779 369 w— 4/3 a7 


é eee 0-876 2 w-2 + O(u-8!8) (2) 
lu 
and 


* ie 
Se = 0806853 w + 350716 wt? — er + 0:00435 w-1+ O(u-4!8) . (3) 


1) University of Kentucky. 
2) Numbers in brackets refer to References, page 63. 
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Table 1 
TE Ena eee” 
lu | a | a* ig — &* - | ; > | < (2, — P*) 
0 | 2-0 3-141 59 
1 2-00862 | 1:90317 | 0:10545 | 333164 | 4:31836 — 0:98672 
2 2.01989 | 1:98755 | 0-03234 | 381493 | 397593 | —90-16100 
3 2.02330 | 201044 | 0-01286 | 4-42973 | 4.40616 | +0-02357 
4 2:02296 | 2:01664 | 0-00632 | 5-10184 | 5-01428 0-087 56 
5 2-02150 | 2.01796 | 0-00354 | 5-80621 | 5-69272 0-11349 
6 | 2:01983 | 2:01765 | 0:00218 | 6:53061 | 6-40692 0.12369 
7 2-01813 | 2:01678 | 0-00135 | 7-26901 | 7-:14223 0-12678 
8 2-:01667 | 2:01577 | 0-00090 | 8-01814 | 7-89125 0-12689 
9 | 2-01536 | 2:01474 | 0-00062 | 8-77507 | 8-64984 0:125:23 
10 2:01428 | 2:01377 | 0-00051 | 9-53817 | 9-41542 0:12275 


It can be seen that in view of the relatively low order of approximation of the 
asymptotic formulae, for « as small as 10 the agreement between asymptotic and 
computed values is surprisingly good. 

Supplementary computations of solutions in the neighborhood of the limit 
cycle solutions were also carried out to investigate numerically the behavior of the 
‘period’ as a function of ‘amplitude’. When the true amplitude a, had been deter- 
mined, additional solutions were obtained for the initial conditions ¥ = 0 and 
Xq = 4 = 0:9 ay, 0-95 ay, 1:05 ay, and 1-1 a, The first subsequent zero of 7 was taken 
to define the ‘semi-period’ to be associated with the ‘amplitude’ @ just as the true 
semi-period P,/2 is associated with the true amplitude a). Table 2 summarizes these 
results. 

The content of Table 2 may be brought out more clearly by using it to construct 
the quantity y = P/P, — 1 as a function of z= 20 (a/a, — 1). The data are sufficient 
that one may fit to them a quartic in z. Thus, let 


y=arz+ P24 yz + 623, 


y being zero at z= 0. The coefficients «, 8, y, and 5 computed from Table 2 are 
given in Table 3 as a function of u. Because « is positive (and in all cases larger than 
}), it is evident that larger initial amplitudes than the limit cycle amplitude give 
longer periods. Moreover, because « increases with mw, a given increase of initial 
amplitude produces a greater increase in period for large w than for small. 


Ee eee 


This result | implies consequences of interest for mechanical systems whose 
motion is described by the van der Pol equation. Thus, suppose that such a system — 


is subject to occasional impulsive perturbations. The effect of such perturbations 
is to cause the system to describe a solution associated with some new initial 
amplitude; the system moves abruptly from one solution of the family of solutions 
whose parameter is a, to some other solution whose parameter is ay + Aa. The first 
point to be noted is that even if the frequency distribution in Aa is symmetric 
because f > 0 there is in the long run a greater lengthening of the period than 
shortening. The mean period of the perturbed system is therefore greater than the 
‘true’ period as given by the appropriate value of LL. 
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Table 2 
P/P, 
aay 
0-90 0-95 1-05 110 
7) 
1 0-928 842 0-962 222 1:042 216 1:088 842 
2 0-857 791 0-925179 1-081 959 1:170813 
3 0-811 063 0-900 063 1-107 907 1:224469 
4 0-779 769 0-884685 1:125451 1-260 820 
5 0-757 797 0-873 236 1-137 823 1:286481 
6 0-741 754 0-864 865 1:146879 1-305 274 
7 0-729642 0-858 546 L535 727 1-319491 
8 0-720 264 0-853 645 1:159044 1:330534 
9 0-712819 0-849755 1-163 264 1:339299 
10 0-706 802 0-846 609 1:166678 1-346 393 
Table 3 
le | 10% | 108 B 108 y | 106 6 


1 | + 39996 
2 78435 3567 —45 42 
3 103513 4432 — 40 2 
4 120423 5066 — 40 2 
5 132334 5528 = 2 
6 141049 5870 my Z 
¥ 147633 6135 — 43 ay: 
8 152743 6343 —44 2 
9 156799 6508 —45 2 

10 | +159989 z 


Secondly, the increase of f with yw implies that this effect is exagerated with 
increasing u. This has bearing, for example, upon systems such as the sequence of 
unstable stars, ranging from the very regular short-period cepheids at one extreme 
to the very irregular and extremely long-period novae at the other. If we assume 
that their radial oscillations are described to a first approximation by the van der 
Pol equation we may conclude that though small « corresponds to the short- 
period cepheids and large mw to the novae, the parameter value associated with the 
latter may be considerably less than is required by the value of the observed mean 
period alone. 

Thirdly, the ‘regularity’ of any such mechanical systems must be greater when 
u is small than when w is large since a given Aa has less effect on account of the 
smaller «. This is reflected in nearly perfectly cyclical behavior of cepheids, the 
approach to the irregularity of relaxation oscillations by the U Geminorum stars 


and repeating novae. 
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Figure 1 


Periodic solutions of the Van der Pol equation over half periods. Each curve is identified by ~ 
the respective parameter value. 
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Figure 2 
Derivative of the periodic solutions of 10 
the Van der Pol equation over half periods. 
Curves are identified by their respective 1 
Van der Pol parameters. 
12 
13 
Figure 3 
Lower half of the phase plane diagrams of the limit cycle : 


solutions of the Van der Pol equation for parameter values 
(Tey A Kay 
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Zusammenfassung 


Die Amplitude und Periode der periodischen Losungen der van der Polschen 
Gleichung sind fiir sehr kleine und fiir sehr grosse Werte des Parameters w durch 
asymptotische Entwicklungen annadhernd bekannt. Fir endliche Werte von wu 
werden die Resultate von numerischen Integrationsmethoden mitgeteilt, und der 
Anschluss an die asymptotischen Werte wird diskutiert. Ausserdem werden Sté- 
rungen der Differentialgleichung betrachtet. 


(Received: July 2, 1959.) 


Plate Thermometer 


By Davip Mexsyn, Lendon, Great Britain’) 


1. Introduction 


The problem of a plate thermometer was first considered by PoHLHAUSEN [1]*). 
The equation of energy is reduced to an ordinary linear differential equation with 
known coefficients, and the temperature function is expressed as a double integral. 
This integral was first evaluated numerically by PoHLHAUSEN for several moderate 
values of Prandtl number o, and later by Eckert and Drewitz [2] for many 
values of o, up to o = 1000. 

For moderate values of o (< 100) the numerical evaluations lead to satis- 
factory results, but for larger o the errors begin to increase, and for o = 1000 the 
‘result is grossly in error (up to 20%). The double integral becomes very slowly 


1) Mathematics Department, Royal College of Science. 
2) Numbers in brackets refer to References, page 68. 
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convergent, as o increases, and it would require several hundred terms to give a 
moderately satisfactory result for o = 1000. 

The aim of the present paper is to indicate a simple method of evaluation of 
this integral, which leads to satisfactory results, as follows. 

The independent variable is changed, and the integral is expressed as a power 
series, depending on a; the first few terms (~ 10) are evaluated in the usual way; : 
the remaining terms are transformed by making use of StiRLING’s formula, and — 
the sum is then transformed into an integral by making use of EULER-MACLAUREN’S — 
equation. The final result, which consists of few terms, is expressed in gamma . 
functions. 


2. The Equation of Energy 


Fluid flows past a semi-infinite plate, placed in the direction of the main flow; 
the velocity and the temperature of the main flow are U, and T, respectively; there 
is no heat transfer between the plate and the fluid. 

The problem consists in finding the temperature of the plate, in this case we 
have to take into account the heat generated by friction. 

The equation of energy is [3] 


a0 d0 0?f\2 
dn? +oaf adn + 26 io =0, (2.1) 
where 
Uz Uy \1/2 
Te ~—_- = ° 
Fee iaigco% rth heats) , 
the boundary conditions are ' 
db at 
‘dn a n= 0., 00 at 1 (2.2) 
/(y) is a solution of BLasius’ equation 
[af faa Oy, (2.3) 
and it satisfies the boundary conditions 
M0) = f'(0) = 9,. f'(co) = 1. (2.4) 


The solution of (2.1), which satisfies the boundary conditions, is 


. 7 : 
O(n) = 20] ("" dn f N?-? dn, (2.5) 
Now [4] pit : | 
f= ae Fo), (2.6) 
where 
4 
F(n) = / f(n) an; 
whence 2 
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where 
an a 7 t1-a3.1° 375 12 
Ein) = —— us [. tinal I a* 4 fees 
ou 6! 9! Ay ; 
a= 0:4696. 


3. Transformation of the Integral 
We have to find 6(0), where 


Ce » 


6(0) = 2ca? / a an | gee dn. (3.1) 


A glance at (3.1) shows that the integral becomes very slowly convergent for 
large o. 
To evaluate (3.1), it is transformed by putting 
F(n) =7; (3.2) 


inverting (3.2) and transforming 6(0) to the rt variable, we find 


co Tt 
a oP on AG o—2)r 
G{O)) == 2 oat | a? a ar |e ¥ a dt; (353) 
0 


0 


from (2.7), (3.2) follows [4] that 


(3.4) 


6 1/3 { nis 71/3 77/3 23 10/3 
Bee (;) lt BOW 260 s9100g.- SI 


Making use of (3.4), where only the first two terms are retained, expanding the 
inner integral in powers of t and integrating in gamma functions, we find 


> r (n+) 


n—0 I'(n +1) r3) (3. +1) 


——— 


6(0) = C otl8 


where 


ae Lo eer 
—— — — = 1-088. 
- Pits ie ne) 


The first ten terms can be easily evaluated; to find the sum of the remaining terms 
we transform the coefficients by making use of the asymptotic expansion of the 


zamma function, and then transform the sum into an integral by making use of 
EULER-MACLAUREN’S formula. 


ZAMP XI/5 


— 
~ 
- 
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Consider the coefficient 
P(n + 3) 
3 . 
ayn, = ’ 


: r(5) Piece Sas.) 


(3.6) 


by making use of the asymptotic expansion [5] 
a il 

=: = —1/2 f = gel 3-7 

Gy = V2 x e-ae* eer oars pr (3.7) 


where only the dominant term is retained, we obtain 


1/6 
e1ls +5.)" 
(3.8) 


An : 
2 ; 1 \n+1/2 1 
4/3 pest a ae 
PPE Ts) 2 (Psa) | Or eel 
Now, 
2 \n 2°90) eae 2 
ated | as N33 
(1+ in) exp {mlog(1 + =} e ante (3.9) 
Transforming similarly the remaining factors we find 
oe = 
bi 2 —}; (3.10 
3 r(5) n4l8 


the expression in (3.5) depending on o becomes 
YP 2 hd <= 
(2= )’ ie nee Bis: = (= = i+ at aS . (Sia 
15 (n— 3) (2) ms ies 


Combining (3.10), (3.11) we obtain for the sum (3.5) 


it 2n 4 (Oe ANY 
= het | ( ) Boe 3.12 
sr(2) 15 (o — 2) on oy - ( 
For large o we can put 
o— 2\n 2 
( = ) = exp [ log 1 = mal = Pag ; (3.13) 
where 


2 
t= 
o 


We transform now the sum (3.12) into an integral by making use of EuLER-Mac 
LAUREN’S formula [5]. 


a+rao 
a 


[ Fe) a= 0 [5 Fla) + Ft 0) + F420) 4-45 Flat eal); (3.14 
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in our case m = 1; whence 


i 2n 4 o— 2\s 1 2n 4 
tp wast bt Pa poate — F ee ; 
a a | * 15 (6 — 2) ral o dl eas CAC) a 
ay (3.15) 
o— 2\"% Pri? 4 2 
¢ 2 aioe 4/8 __ yy —7/8 — —1/3| p-«! 
x ( - ) ai (x > @ 2 a Fg we je nan . 
ny 
Integrating by parts, we find 
fo, 2) 
3 3 1 € 
20 ere i a ~ 4+ 5) ee 
a 
! ny 9 ny, 2) (3.16) 
beg € oO 2 2 
ni Hits Baca 2) y (zs em 
where y is the incomplete gamma function 
2 2 
Ps =) = 73) 
finally . 
2 
pao ers) | 
6(0) = 1-088 o1/8/1 +S" 4 
n-1 (3m +1) r(5) (n+ 1) 
o—2\n (2n—3) » Scie 4 a 
(".) a3 fiat : Ss 
(m5) 
1 ae (m ¢ 4 | 3.17 
sigs (3.17) 


2 15 ( A 5) os, 
00 
a + 0-2462 | ; 
ny 
where (3.12) 
0:2462 = 


RSQ t SSNPS PO 
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In the table are given a few of the results obtained (the current values are indicated 
in parentheses) : 


: 
to} | 15 | 100 | 300 | 1000 | co : 


6(0) 3-6 7:8 11-6 17:9 1:9 o/8 
| (7) (10) (14) 
REFERENCES : 
[1] E. Ponrnausen, Z. angew. Math. Mech. 7, 115 (1921). 
[2] E. Eckert and O. Drewitz, Luftfahrtforschung /9, 187 (1942). 
[3] H. Scuricutine, Grenzschicht-Theorie, 3rd ed. (Verlag G. Braun, Karlsruhe 1958). 
[4] D. Mexsyn, Proc. Roy. Soc. [A] 237, 543 (1956). 
[5] E. T. WarttaKerR and G. N. Watson, A Course of Modern Analysis, 4th ed. 
(Cambridge University Press, 1949). 


Zusammenfassung 


Das Doppelintegral, das mit der Temperatur eines Plattenthermometers zu- 
sammenhangt, konvergiert langsam bei grossen Werten der Prandtl-Zahl o. E 
wird eine Methode angedeutet, um dieses Integral fiir alle nicht zu kleinen Wert 
von o auszurechnen. Bei o = 1000 ist der ausgerechnete Wert ungefahr 20% grésser 
als gewohnlich angenommen wird. 


(Received: July 20, 1959.) 
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Réunion d’automne de la Société Suisse de Physique 
dans le cadre de la 139° assemblée annuelle de la Société Helvétiqu 
des Sciences naturelles, 4 Lausanne le 12 septembre 1959 


Compte rendu des communications de physique et de mathématiques appliquées 


Ein stabiler Thermostat zur Messung des Rauschens und des Drift 
von Halbleiter-Schaltelementen. Von E. BALpDINGER und A. Mater, Basel! 
Messungen des Rauschens und des Drifts von Halbleiter-Schaltelementen sin 
sowohl von praktischem als auch von theoretischem Interesse. Praktisch geben si 
Aufschluss iiber die in Schaltungen erreichbare Genauigkeit, theoretisch liefern si 
Aussagen tiber die Halbleiter selbst. Speziell mit den Vorgangen an deren Obe 


1) Institut fiir angewandte Physik, Universitat Basel. 
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flache ist der «excess noise» verkniipft, dessen Spektrum im bisher untersuchten 
Frequenzgebiet (oberhalb etwa 10-4 Hz)-mit (1/f)" verlauft, wobei n, je nach dem 
Objekt, in der Regel zwischen 1 und 1,5 liegt®). Bei geniigend tiefen Frequenzen 
kann man eine Abweichung von diesem Verlauf erwarten3), deren Lage und Art 
weiteren Aufschluss iiber das Oberflachenrauschen liefern wiirde. 

Entsprechende Messungen von Rauschspannungen erfordern, unter Beriick- 
sichtigung des Temperaturkoeffizienten, dass die Temperatur der Objekte iiber 
mehrere Tage auf etwa 1mGrad konstant gehalten wird. Im folgenden soll ein 
Thermostat beschrieben werden, der diesen Anforderungen geniigt und der auch 
ganz allgemein als Prazisionsthermostat verwendbar ist. 


Aufbau 


Der thermostierte Raum besteht aus einem geriihrten 20-Liter-Olbad, dem 
innern Thermostaten, der warmeisoliert (0,5 W/Grad) in einem zweiten, dem 4us- 
sern Thermostaten, eingeschlossen ist. Letzterer reduziert die Schwankungen der 
Raumtemperatur um den Faktor 400, wird aber durch elektrische Zuleitungen und 
die mechanische Antriebswelle (Plexiglas) zwischen Aussenraum und innerem 
Thermostaten zu einem geringen Teil tiberbriickt. 

Die Festlegung der Temperatur des inneren Thermostaten erfordert die Mittel 
der Prazisionsthermometrie. Abgesehen von den Temperaturfixpunkten, die sich 
fiir Thermostatregelung nicht unmittelbar eignen, gelten als genaueste Instrumente 
im fraglichen Temperaturgebiet um 30°C die Platin-Widerstandsthermometer. 

Im vorliegenden Thermostaten wird ein 100-Q-Pt-Thermometer verwendet, das 
in einer Briickenschaltung mit 0,1 V Wechselspannung von 1000 Hz gespiesen 
wird. Die Ausgangsspannung der Briicke wird verstarkt, in einer dem Ringmodu- 
lator ahnlichen Schaltung gleichgerichtet und zur Steuerung der Heizung beniitzt. 
Die gesamte Schaltung hat eine Bandbreite von etwa 0,5 Hz. Sie ist voll transisto- 
risiert. 


Temperaturkonstanz 


Die verstarkte und gleichgerichtete Ausgangsspannung der Briicke ist ein 
Mass fiir die Temperatur des Thermostaten. Aus ihrer Registrierung lasst sich 
dessen zeitliche Konstanz ermitteln (Tabelle 1). Allerdings bleiben auf diese Art — 
gewisse Effekte, wie etwa eine Veranderung des Pt-Thermometers selbst, unbemerkt. 
Solche Erscheinungen kénnen nur durch Verwendung mehrerer voneinander unab- 
hangiger Temperaturregistrierungen aufgedeckt werden. 


Tabelle 1 
Tempervaturschwankungen im Thermostaten 


Quadratischer Mittelwert Maximale Abweichung von 
ena eiaest suet der Temperaturschwankung | der mittleren Temperatur 
10 min 0,02 mGrad -- 0,06 mGrad | 
100 h Dauerbetrieb 0,06 mGrad + 0,15 mGrad 


R. H. Krvcsron, J. appl. Phys. 27, 114 (1956). : 
3) A. L. McWuorter, Semiconductor Surface Physics, edited by R. H, Kryesron (Philadelphia 
S 
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Das Temperaturfeld im innern Thermostaten ist notwendig inhomogen, sobald 
in dessen Innenraum Warme zugefiihrt wird. Bei einer Warmezufuhr von 1,8 W 
und einer Riihrgeschwindigkeit von 130 U./min-ergab sich in einer provisorischen 
Messung eine Temperaturdifferenz von 3 mGrad zwischen Zentrum und Rand. Bei 
der Untersuchung von Objekten kleiner Leistung wird man normalerweise mit 
etwa einem Zehntel dieser Heizleistung auskommen und infolgedessen nur mit 
0,3 mGrad Temperaturdifferenz rechnen missen. 

Ein stdérender Nebeneffekt des Riihrens ist zu beachten. Der Warmeiibergang 
von den geheizten Elementen an das Olbad haingt von der Grésse der Oberflache 
und der Stromungsgeschwindigkeit des Ols ab. Letztere ist durch das Rihren 
Schwankungen unterworfen, welche sich, besonders stark bei kleiner Oberflache, in 
Schwankungen des Warmeiibergangs und, damit der Temperatur der geheizten 
Elemente auswirken. So resultiert zum Beispiel in einer Intermetall-Zener-Diode Z6 
bei 60 mW Leistung eine mittlere Temperaturschwankung von 5 mGrad. Dieser 
Effekt lasst sich dadurch stark reduzieren, dass die Zener-Diode mit Woodschem 
Metall in einen 20-cm’-Messingzylinder eingegossen wird. Man profitiert dabei sowohl 
von der grossen Oberflache als auch von der 6rtlichen und zeitlichen Mittelung. 


Experimente 


Im Rahmen unserer Untersuchung iiber stabilisierte Transistornetzgerate wurden 
die temperaturempfindlichen Teile einer Querschaltung in den beschriebenen Ther- 
mostaten gebracht. Die Schaltung ist eine fiir diesen Versuch leicht modifizierte 
Ausfihrung (Figur 1) eines bereits beschriebenen*) Gerates. Die Abanderung besteht 


Speisung 
12V stab. 


Speisung 


77V stab. 


Querschaltung 
auf £07 °C thermostiert 


Referenz- 
Spannung 


Z5 


Symm. 
Verstarker 


2x00 76 


im beschriebenen Thermostaten 


Ld 


Figur 1 
Blockschema des untersuchten Netzgerates. 


Ausgang 
6V/2A 


in der Hauptsache darin, dass die Referenzspannung direkt von einer Zener-Diode 
Z 5 abgenommen wird, die ihrerseits iiber eine Zener-Dioden-Vorstabilisierun aus 
einer hinreichend stabilisierten Spannungsquelle gespiesen wird. — - 
sind in Tabelle 2 zusammengefasst. ' 


4) E, BALDINGER und W, Czaja, Z. angew. Math. Phys. 9a, 291 (1958). 


Die Ergebnisse 


: 
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Tabelle 2 


71 


Rauschen und Drift des stabilisierten Tvransistornetzgerdtes nach Figur 1, sowie des einge- 
bauten symmetrischen Gleichstromverstdvkers. Zum Vergleich sind entsprechende Werte fiir 


einen Autoakkumulator angegeben 


Gleichstromverstarker 


Akkumulator 


N ceratausgane 
ae Nee allein, auf dessen 6 V / 80 Ah, 
Eingang bezogen unbelastet 
Rauschspannung 
quadratischer Mittelwert 
liber 650s mit 1 Hz 
oberer Bandgrenze . 2 LV U2 NG <0,3 wV 
Drift tiber 1 Monat. 1,6 uV/h 0,6 uV/h 150 wV/h 
n=7,65 +1000 
n=155 
Faktor 40 aed 
lin 


He tees 7 tH 
n=100 i 7 ~ 
ps 3 
L aq i — 1 
1000 100 10 1 gth 


Figur 2 
Rauschmessungen bei extrem tiefen Frequenzen: b Zener-Diode Z5 der Referenzspannung; 
PY Transistorpaar OC 76 aus dem symmetrischen Verstarker (Messungen auf dessen Eingang 
bezogen). U ist die quadratisch gemittelte Rauschspannung, gemessen wahrend der Zeit T mit 
einer oberen Bandgrenze von 1 Hz. Die Bezugsspannung U, berechnet sich aus dem «excess noise» 
bei 470 Hz (vel. Figur 3). Die ausgezogenen Linien geben, mit verschiedenem n, den Verlaut von 
(U/U,) fir ein (1/f)"-Spektrum, das bei 470 Hz an die Tonfrequenzmessungen anschliesst. Die 
angegebenen Fehlergrenzen sind mittlere Fehler. 
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Das Rausch-Spektrum von Halbleiter-Schaltelementen 


Die erwahnten Experimente waren so angelegt, dass sie direkt Aussagen iiber das 
Rauschen der Zener-Diode Z 5 (Referenzspannung) und der beiden Transistoren 
OC 76 (Eingangsstufe symmetrischer Verstarker) bei tiefen Frequenzen zuliessen 
(Zener-Strom 5 mA; Emitterstrome je 0,6 mA). Diese Aussagen wurden mit 
Rauschmessungen im Tonfrequenzgebiet .verglichen (Figuren 2 und 3). 


10 -15 |_ 


—t SS 


10° 10° 10° Hz 


f —> 
Figur 3 
Rauschspektren im Tonfrequenzgebiet: © Zener-Diode Z 5 der Referenzspannung (einschliesslich 
weisses Rauschen), Uy = 0,04 wV; : Transistorpaar OC 76 aus dem symmetrischen Verstarker 
(Messungen auf dessen Eingang bezogen, weisses Rauschen subtrahiert), Uy) = 0,59 UV, m = 1,65; 
Q weiterer Transistor OC 76 (weisses Rauschen subtrahiert), bei 470 Hz auf gleiche Rauschspan- 
nung normiert wie oe 
Die Ordinate dU?/df ist die spektrale Dichte des Quadrates der Rausch-Spannung. 


dU2 
U, = 
M (\/ df ADeS 


charakterisiert die Grésse des «excess noise» bei 470 Hz. Die angegebenen Fehlergrenzen sind 
mittlere Fehler. 


Zur Auswertung der Messungen bei tiefen Frequenzen ist folgendes zu bemerken. 
Es wurde jeweils wahrend eines Zeitraumes T die entsprechende Spannung wieder- 
holt in etwa gleichen Zeitabstanden gemessen und durch einen geeigneten linearen 
Gang approximiert. Die Abweichung der Messpunkte von der resultierenden Gera- 
den wurde daraufhin quadratisch gemittelt. Die Wahl dieser Auswerteart geschah 
in der Meinung, dass dabei sowohl ein Drift als auch die im Intervall T vorwiegend 
linearen Beitrage von Frequenzen < 1/2 T weitgehend verschwinden, so dass ma 
das so ermittelte Rauschen mit einer fiir vorliegende Zwecke hinreichenden Ge 


Vol. XI, 1960 Varia — Miscellaneous — Divers . 73 


nauigkeit dem Integral iiber das Rauschspektrum von hp yan biss, = 1rlz 
gleichsetzen darf. 

Die Messungen an der Zener-Diode deuten auf die Existenz eines einheitlichen 
(1/f)"-Spektrums mit = 1,55. Die Messungen an den Transistoren zeigen nicht 
dasselbe einheitliche Bild. Sowohl im Tonfrequenzgebiet als auch bei extrem tiefen 
Frequenzen verlauft das Spektrum wie (1/f)" mit m = 1,65, doch differieren die 
absoluten Betrage der entsprechenden Rauschquadrate um den Faktor (40)? = 1600. 
Der Ubergang, der nach diesen Messungen zwischen 0,5 kHz und etwa 1 mHz 
stattfinden muss, ist in der Messung eines in der Schaltung nicht verwendeten 
OC 76 bei 250 Hz und 75 Hz angedeutet. 

Fir die Rauschmessungen im Tonfrequenzgebiet sei Herrn BitGerR unser herz- 
licher Dank ausgesprochen. Der « Kommission zur Férderung der wissenschaftlichen 
Forschung» des Bundes danken wir fiir die Finanzierung der Untersuchungen. 


Uber die Méglichkeit digitaler hydraulischer Elemente. Von H. GLArtit 
und H. R. MUxier, Ziirich). 

Elementare Uberlegungen zeigen, dass hydraulische Schaltelemente Ansprech- 
zeiten der Gréssenordnung Millisekunden gestatten. Hierbei sind ein Druck von etwa 
5 atii, ein Hub von 2mm und eine Lange der beschleunigten Saule von 20 cm an- 
genommen. 

Figur 1 zeigt 2 hydraulische Multivibratoren, die sich nur beziiglich Grésse, 
nicht aber durch ihr Arbeitsprinzip unterscheiden. Die Frequenz des kleinen 
Multivibrators liegt etwa 4mal hoher als die des groésseren. 


Figur 1 
Geometrisch ahnliche hydraulische Multivibratoren. 


Die Frequenz eines Multivibrators stellt sich so ein, dass die Verstarkung 1 
betragt. Als Verstarkung bezeichnen wir das Verhaltnis zwischen aufzuwendender 
und steuerbarer Energie in einem bestimmten Zeitintervall. In einem Fall betrach- 
ten wir die Energie, die notwendig ist, um ein Ventil innerhalb einer vorgegebenen 
Zeit (Ansprechzeit) auszusteuern. Im zweiten Fall betrachten wir die Energie, die 
von einem ausgesteuerten Ventil in derselben Zeit kontrolliert werden kann. 


5) IBM-Forschungslaboratorium. 
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Ventil Relais inacgein 
Wahrheitstabelle 
eau M 
; So tI} No——+ 
Oy etal as B 
y Y Hl ate ee B : 
Lire: A 
X=M(AB+A'B') 
X=M(A'+B) + N(AB') N(A'B+AB') X=A(B+C) 

Figur 2 


Hydraulisches logisches Element und entsprechende elektrische logische Elemente. 


Figur 2 zeigt ein hydraulisches logisches Grundelement und die entsprechenden 
elektrischen logischen Grundelemente Transistor und Relais. Die Wahrheitstabelle 
definiert die Stellung des Ventils als Funktion der Drucke A und B. 

Durch Fixieren oder Parallelschalten von einzelnen Eingangen kénnen Unter- 
funktionen mit (n — 1), (n — 2), ... usw. Variabeln erzeugt werden. Diese Methode 
gestattet mit ein und demselben Element mehrere logische Funktionen zu erzeugen, 

Die folgende Tabelle zeigt die Anzahl Unterfunktionen verschiedener Elemen- 
tarschaltungen. 


| Transistor Ventil Relais 
Funktionen mit 4 Variablen — 1 1 
Funktionen mit 3 Variablen 1 (4) ie 7 
Funktionen mit 2 Variablen 3 (4) 4 6 


Die eingeklammerten Gréssen in der Transistorkolonne haben nur Bedeutung, 
wenn wir wahlweise das eine oder andere Grundelement zulassen (pnp- oder npn- 
Transistor, Kollektor- oder Emitterbasisschaltung). Hingegen kommt ihnen nur 
wenig Bedeutung zu, wenn wir Mikroprogrammierung vorsehen. Der Transistor ist 
dem hydraulischen logischen Element und dem Relais beziiglich Schaltzeit einige 
Grossenordnungen iiberlegen. Der Anwendungsbereich hydraulischer logischer 
Elemente diirfte sich deshalb auf Anwendungen mit relativ geringer Steuerge- 
schwindigkeit beschranken. 

Die Frage der Miniaturisation betrifft vor allem die Schaltzeit, das Volumen 
und die Antriebsenergie. Die Reduktion des Volumens ist offensichtlich, wahrend 
Antriebsenergie und Schaltzeit wesentlich von den Voraussetzungen, wie Viskositat 
und Toleranzen, abhangen. Die Untersuchungen beschranken sich auf die Schalt- 
zeit. Den Berechnungen wurde die Annahme zugrunde gelegt, dass alle 3 Dimen- 
sionen um den gleichen Faktor « reduziert werden. Vernachlissigen wir die Reibung 
so ist die Zeitreduktion proportional «, wie eine einfache Rechnung zeigt. Fiihren 
wir die Reibungen ein, so ist das Problem nur noch numerisch lésbar. Die nume- 
rischen Rechnungen wurden auf einer IBM 650 durchgefiihrt. 


Figur 3 zeigt den Zeitreduktionsfaktor y (Verhaltnis zwischen neuer Schaltzeit ¢ 
und urspriinglicher Schaltzeit ¢)) in Abhangigkeit von «. 
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Der Einfluss der Reibung macht sich erst bei der fiir Hydraulikanwendung 
ungewohnlich hohen Zahigkeit von 100°cP bemerkbar. Bemerkenswert ist dabei 
noch die Tatsache, dass die Reibungsverluste in der Scherstromung zwischen Kol- 
ben und Bohrung die viskosen Verluste der Laminarstromung um mindestens eine 
Groéssenordnung iibertreffen. 

Zusammentassend kann folgendes gesagt werden: 

Hydraulische logische Elemente gestatten Schaltzeiten der Gréssenordnung 
Millisekunden. Sowohl logisch wie auch beziiglich der Betatigungsenergie besteht 
grosse Ahnlichkeit mit Relais. Miniaturisation erhdht die Schaltgeschwindigkeit. 
Die Verwendbarkeit hoch viskoser Fliissigkeiten verringert weitgehend die auftre- 
tenden Toleranzprobleme. 


? 


net und: 10 cP 7:1000 cP 


Figur 3 


Zeitreduktionsfaktor y als Funktion von «. Druck | ati, Parameter: Hub. 1 6,4mm; 2 3,2 mm; 
31,6mm; 4.0,8mm; 4 0,4 mm. 


Détermination des paramétres h,, et h., de transistors. Par R. Dessov- 
LAvy et B. SEcRETAN, Lausanne’). 


I. Introduction 


er 


1. Symboles et schéma équivalent du transistor a jonction dans le montage’ 
émetteur commun (figure 1) : 
2. Définitions: 


En trés basse fréquence (~ 100 Hz) ces 2 grandeurs sont réelles et valent 


a roximativement: 
i = IR fe 1 S Ry 
hy, = > fax Fe ie 


6) Ecole Polytechnique de l'Université de Lausanne. Laboratoire d’Electronique. 


es 
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Nous montrons dans ce qui suit la possibilité de déterminer ces valeurs par des 
mesures A une fréquence supérieure (~ 1 kHz) grace a une méthode de neutralisa- 
eee Méthodes de mesure. Nous ne donnerons sur les figures 2, 3 et 4 que le circuit 
équivalent pour l’alternatif. Il va de soi que pour les mesures, des sources de 
polarisation adéquates doivent étre prévues. oa ay ; 

a) Mesure directe des courants et tensions alternatives I,, U, et U, selon la 
figure 2 (pour la mesure I, Ry = oo, pour celle de U, Ry = 0). 

Les courbes a des figures 5 resp. 6 montrent la dépendance de /y, resp. Axe 
fonction de la fréquence. On constate qu’a 1 kHz, méme déja a 300 Hz, on mesure 
des valeurs nettement différentes de celles obtenues en trés basse fréquence 
(~ 100 Hz). Ceci s’explique par le couplage capacitif entre B’ et C (cf. figure 1) 
di a Cy. 


Figure 1 Z, 


Transistor 4 jonction: définition des courants et tensions et schéma équivalent pour de faibles 
signaux. 


b) Mesure en pont a 1 kHz selon les figures 3 et 4. 

Les courbes a des figures 7 resp. 8 représentent la tension de deséquilibre Up 
du pont en fenction des résistances d’équilibrage R, resp. R,. On constate qu’a la 
fréquence utilisée de 1 kHz la tension Up ne s’annule pas mais passe par un mini- 
mum d’ailleurs trés peu marqué. Ceci provient du fait que U, et U, ne sont pas en 
phase dans la mesure de h,, et que l’admittance inconnue hy est complexe. 


II. Neutralisation 


1. Principe: La neutralisation consiste a faire circuler dans R, et Z, un courant 
de phase opposée a celle du courant dfi A Cy. La méthode classique utilisée ici est 
de connecter une capacité de neutralisation C,, (C;,) entre la base et une source de 
tension — U, (en pointillé sur les figures 2, 3 et 4) 

2. Détermination de h,,. Hypothéses: 


= il 
a) 45<2Z 5; b) Ry te Zalk peer, C)RouC ene alle 


n 
On trouve pour hy: 
> R, 
hie Spy LL + OF R, Rs Cy (Cz — Cy) + jo [RC, —C, (Ry + Ry) — oat) ; 


3 
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Figure 2 


Mesure directe des paramétres hy) et Aig9. En pointillé: neutralisation. 


Figure 3 


Détermination de h,, par une mesure en pont. En pointillé: neutralisation. 


Figure 4 


Détermination de /,, par une mesure en pont. En pointillé: neutralisation. 


hab 
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La neutralisation aura lieu lorsque 


= R C, Ry 
Im {hy} =9, Cnr= "Fok Cs (1 ‘ ee k 


eer 


Pour cette valeur de C,, 


r 
| 


a condition que 
é= w? R, BR, C, (C, —C,) <1. 


Sinn 


3. Détermination de hyy. On a facilement j 


2 1 oe 
hos = p- +S te 


4 U, 


| 
= é 
En admettant les mémes hypothéses que pour /,, on trouve: 
: 
. 
. 
- 


: 


a condition que 


v=o? RRC, (Cre OC) <i. 


2 1 R, ise RS 
Io & —-t+S>*[1 +o FR, RC, (C3; —C,)] +7 @S [(cs - Ca) Rg CR - 
Ry Tes Tis 
4 
La neutralisation aura leu lorsque : 
Lut {h.9) = Oo 
dot | 
Ci k 
c.=C (: ee | 
n 3 Cr R; 
a leur pour laquelle h,, vaudra : 
1 Hie 
Noo => WER aii S wie 


III. Résultats expérimentaux _ 


Les mesures ont été effectuées sur des transistors OC 71 au germanium du es 
courant. Pour 7,=1mA et u,,= —5 V on a obtenu avec un exemplaire plutot 
défavorable les résultats suivants: 


Sur la base de mesures complémentaires nous avons déterminé pour ce tran- 
sistor les éléments de la figure 1: 


Notas hee dad Rg. = 558 MQ, R,= 200k09 
ye ol aie g = 69 pF, C, = 56pF, Ce. a5 
S =40mAJV. . ae 
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6 
10° h,, 


Figure 5 


Variation de fy, en fonction de 
la fréquence. Mesure selon 
schéma de la figure 2. @ sans 
neutralisation; b avec neutrali- 

sation. 10 20 50 100 200 500 1000 3000 Hz 


palvT hoo 
60 


40 


Figure 6 20 


Variation de hy, en fonction 

de la fréquence. Mesure selon 

schéma de la figure 2. a sans 

neutralisation; bavecneutrali- 0 
sation. 


ri Figure 7 


Mesure de fy, au pont selon 
le schéma de la figure 3. Ten- 
_ sion Up en fonction de la ré- 


2 sistance d’équilibrage 
R, (R, = 100kO; U, = 0,5 V). 
. b a Mise” ? 
a sans neutralisation; b avec 0 20 f 5 40 60 80 100 2 
neutralisation. 10 hie 
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pV 


1000 \ 
\ 
750 \ 
\ 
500 \ 
\ 
“Ss 
250 b 
\ oe 
Ne Ae R 
0 eon c 
0 20 40 60 80 100 kQ 
tng 


Figure 8 
Mesure de fig) au pont selon le schéma de la figure 4. Tension Up en fonction de la résistance d’équi- 
librage R, (Rg = 100 Q; U, = 0,5 V). a sans neutralisation; b avec neutralisation. 


IV. Conclusion 


Les mesures confirment le fait que l’on peut mesurer a 1 kHz grace a la neutrali- 
sation les valeurs recherchées de hy et hy. pour les trés basses fréquences. Les erreurs” 
systématiques « et e’ restent dans des limites encore acceptables (« = 6% ete’ = 5%). 
Enfin cette méthode de neutralisation permet de calculer C, a partir de C, si lon 
connait ht, io,dt. et Cp. 


' 
; 
{ 
: 


International Symposium on Active Networks and Feedback Systems 


New York, April 19-21, 1960 | 


The Microwave Research Institute of the Polytechnic Institute of Brooklyn 
announces that the tenth of its annual International Symposia will be held at the 
Engineering Societies Building, 33 West 39 Street, New York, New York, Tuesday 
through Thursday, April 19-21, 1960. The papers given will include the blow 
topics: 

Physical Realizability of Active Networks, . 
Representation and Synthesis of Active Networks, 

Synthesis of Feedback Networks, : 
Time Varying Systems, 
Applications of Active Network Theory to Electron Devices, | 
Applications of Active Network Theory to Feedback Systems. | 


Correspondence should be addressed to: Professor HERBERT J. CARLIN, Micro- 
wave Research Institute, 55 Johnson Street, Brooklyn 1, New York. 


JEROME Fox 
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Festschrift fir Alexander M. Ostrowski. 1. und 2. Teil. Archiv der Mathe- 
matik 5, Fasc. 1-3 und 4-6 (Verlag Birkhauser, Basel und Stuttgart 1954). 521°S., 
40 Abb.; Fr. 60.-. 

Der fiinfte Band des Archives der Mathematik erschien — einem schénen Brauch 
folgend — als Festschrift zu Ehren von Prof. ALEXANDER M. Ostrowskt, der am 
25. September 1953 seinen 60. Geburtstag beging. Eine grosse Zahl von Freunden, 
Kollegen und Schiilern aus der ganzen Welt widmeten dem verdienten Mathema- 
tiker und geschatzten Lehrer kleinere oder gréssere Abhandlungen aus den ver- 
schiedensten Gebieten der Mathematik und demonstrierten damit auch das um- 
fassende Wirken des Jubilars. Die Arbeiten entstammen vorwiegend der reinen 
Mathematik; doch finden sich darunter auch eine Reihe von Untersuchungen aus 
der angewandten Mathematik oder solche, die direkt oder indirekt zu ihr Bezug 
haben und damit speziell auch den Leser der ZAMP interessieren werden. Die 
beiden Teile dieser ausgezeichneten Festschrift umfassen 70 Arbeiten, so dass es 
unmoglich ist, hier auf sie einzeln einzugehen. E. Rotu-DESMEULES 


Applied Mathematics and Mechanics (Prikladnaia matematika i mekha- 
nika). Translation Editor: GEORGE HERRMANN (Pergamon Press Inc., New York und 
London). Erscheint zweimonatlich, jahrlich etwa 1200 Seiten. Jahresabonnement 
$ 35. ; 
Es ist bekannt, dass die russische wissenschaftliche Literatur — abgesehen von 
den Referaten — in der westlichen Welt vielfach wenig oder verspatet beachtet wird. 
Der Grund fiir diesen Zustand ist wohl in erster Linie in den sprachlichen Schwierig- 
keiten («language barrier») zu suchen; ausserdem sind manche Zeitschriften nur 
schwer zugdnglich. Dabei ist die Zahl der hier interessierenden russischen Zeit- 
schriften, Nachrichtenblatter und Berichte beachtlich: ein halbes Dutzend fiir 
Mathematik, ein Dutzend fiir Physik, etwa drei Dutzend fiir Technik (zum Teil sehr 
stark spezialisiert), wozu noch eine Anzahl von Zeitschriften in einer der Sprachen 
der Unionsrepubliken kommen. Fiir den Forscher und Wissenschaftler im Westen 
‘stellt es jedenfalls eine erhebliche Erleichterung dar, wenn ihm wenigstens die — 
wichtigsten Zeitschriften in einer Ubersetzung zuganglich sind. 

Zu den bedeutenden russischen Zeitschriften geh6rt die « Prikladnaia matematika 
imekhanika», die nun seit 1958 vom Pergamon Institute (New York) im Auftrage 
der American Society of Mechanical Engineers und mit Unterstiitzung der National 
Science Foundation iibersetzt und herausgegeben wird. Diese Zeitschrift, von der 

gegenwiartig (1959) der 23. Band erscheint, umfasst das Gebiet der angewandten 

‘Mathematik und der Mechanik (fester, fliissiger und gasformiger K6rper), entspricht 

also etwa der ZAMP. Um einen Einblick in das Arbeitsgebiet und die behandelten 

Probleme zu geben, seien hier wenigstens die Titel der zum Teil recht bemerkens- 
‘werten Arbeiten der vorliegenden ersten Nummer der Ubersetzung (Band 22, Nr. 1) 
-vermerkt: 

_ N.N. Kocuina and N. S. Mev’nrkova, Strong Point-Blasts in a Compressible 
“Medium. V. A. Ecorov, On the Solution of a Degenerate Variational Problem and the 
Optimum Climb of a Cosmic (Space) Rocket. Ia. KurTSVEIL’, On Generalized Ordinary 
Differential Equations Possessing Discontinuous Solutions. V. 1. ZuBov, On a Method 
‘of Investigating the Stability of a Null-Solution in Doubtful Cases. 1. B. CHELPANov, 
Oscillations of a System with a Relay of Advancing Characteristics. B. A. BERG, 


ZAMP X1/6 


— 


. 
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Deformational Anisotropy. lu. 1. KADASHEVICH and V. V. NovozuiLov, The Theory | 
of Plasticity Which Takes into Account Residual Microstresses. D. D. IvLEv, On the 
General Equations of the Theory of Ideal Plasticity and of. Statics of Granular 
(Pulverulent) Media. V. D. KLIUSHNIKOV, On Plasticity Laws for Work-Hardening 
Materials. L. A. SHapovatov, The Influence of Non-Uniform Heating on the Stability 
of a Compressed Bar. V. 1. Mossakovskul, Pressure of a Circular Die (Punch) on an 
Elastic Half-Space, Whose Modulus of Elasticity Is an Exponential Function of 
Depth. V. V. Ketpysu, Application of Slender Body Theory to the Calculation of 
Aerodynamic Properties of Low Aspect Ratio Wings with Nacelles at Their Tips. 
G. N. Kopyiov, Aerodynamic Characteristics of Thin Wedge-Shaped Profiles im 
Trvansonic Flows. S. Manotov, A Special Case of the Existence of Small Periodic 
Motions of two Penduli, Subjected to Uniform Rotation. G. K. POZHARITSKII, On a@ 
Property of the First Approximation System. { 
Die Durchsicht der Arbeiten der vorliegenden Nummer bestatigt, dass im all- 
gemeinen der russische Wissenschaftler iiber die westlichen Arbeiten eingehend 
informiert ist (eine Beziehung, die nicht ohne weiteres umkehrbar ist, aber auch 

verschiedene Interpretationen zulasst). Das Verhaltnis der zitierten russischen un 
fremdsprachigen Arbeiten ist etwa 9:8. — Die englische Ubersetzung der PMM wird 
jahrlich ungefahr 1200 Seiten umfassen, was den etwas mehr als 850 Seiten des 
russischen Originals entspricht. Ohne Zweifel wird diese Publikation auf ein erheb- 
liches Interesse stossen. Erwahnt sei noch, dass von jeder Arbeit auch Sonderdrucke 
zum einheitlichen Preis von 2.00 $ (10 s) erhiltlich sind, was sicher sehr zu begriissen 
ist in Anbetracht des doch recht hohen Preises eines ganzen Jahrganges. ; 
E. RotH-DESMEULES 


Anwendung von Rechenmaschinen bei der Berechnung von Regely 
vorgangen (Vortrage gehalten bei einer Tagung des Fachausschusses Regelungs- 
mathematik der GAMM in Diisseldorf, 8. November 1957). Herausgeber: W. Op- 
PELT (Oldenbourg, Miinchen 1958). 128 S., 121 Abb.; DM 16.80. ; 

Nachdem sich in den letzten Jahrzehnten die Regeltechnik immer starker zur 
Lésung mathematischer Probleme der programmgesteuerten Digitalrechengerate 
und auch der Analogrechner bediente, widmete die GAMM eine Tagung der Behand- 
lung von Fragen, die sich bei der Anwendung moderner Rechenmaschinen in der 
Regelungsmathematik ergeben. In sieben ausgezeichneten Vortragen kamen dabe 
sowohl Probleme aus der digitalen wie auch aus der analogen Technik zur Sprache, 
welche dann in der Reihe « Beihefte zur Regelungstechnik» des Oldenbourg-Verlages 
verdffentlicht wurden. ; 

E. Bukovics (Wien) gibt einen Uberblick iiber die Einsatzméglichkeiten digi- 
taler Rechengerate zum Lésen von Grundaufgaben der Regelungstechnik. Er zeigt 
wie mit relativ kleinen digitalen Maschinen Stabilitats- und Optimierungsberech: 
nungen von Regelkreisen durchgefiihrt werden kénnen. 

J. B. Reswick (Cambridge, Mass., USA) erlautert eine einfache graphische 
Methode, um aus einer gegebenen Autokorrelationsfunktion der Eingangsgrésse 
eines Regelkreises und der Crosskorrelationsfunktion zwischen Eingang und Aus- 
gang die Ubertragungsfunktion des Regelkreises zu finden. | 

; A. LEONHARD (Stuttgart) gibt einige Anwendungen eines elektronischen Ge 
rates an, welches gestattet, den komplexen Wert von Polynomen zu berechneil 


Besonders bei der analytischen Behandlung von Regelauf ist si i 
Gerat als gutes Mittel zur Zeitersparnis. Sloot ak ee 


Bis heute sind schon einige 
optimal zu programmieren. R. H 


al 


Verfahren bekannt, Analogierechenschaltunge 
ERSCHEL (Ulm) zeigt eine weitere Programmie- 
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rungsart, die sich insbesondere bei der Behandlung nichtlinearer Probleme sehr gut 
bewahren soll. 

Ein wertvoller Vergleich der Berechnung eines Regelvorganges mit Rechen- 
anlagen verschiedener Typen (Digitalrechner IBM 650, elektronischer Analogie- 
rechner, mechanische Integrieranlage) stammt von den Herren O, FOLLINGER und 
G. SCHNEIDER (Frankfurt). 

E,. BUHLER (Mannheim) verfasste einen Beitrag iiber die elektronische Nach- 
bildung eines mechanischen Systems mit Reibung. Sorgfaltig werden besonders die 
Reibungsverhiltnisse beim Ubergang von Haftreibung zu Gleitreibung nachgebildet, 
wobei allerdings der fiir diese Nachbildung betriebene elektronische Aufwand frag- 
wiirdig gross ist. 

Ein letzter Vortrag von H. WiTsENHAUSEN (Briissel) iiber die Anwendung von 
Analogierechenmaschinen auf die Optimierung von unstetigen Regelkreisen mit 
statistisch schwankenden Eingangsgréssen behandelt ein heute sehr aktuelles 
Thema, da man in jiingster Zeit immer mehr danach trachtet, Regelkreise nicht fiir 
Eingangsgréssen, die als Sprungfunktion auftreten, zu optimieren, sondern fiir 
Statistisch schwankende Eingangsgréssen, da die wirklich vorkommenden Stér- 
funktionen meist auch nur statistisch erfasst werden kénnen. 

Einige interessante Diskussionsbeitrage bilden den Abschluss dieses kleinen 
Werkes, das jedem, der sich mit Regeltechnik und Analogierechentechnik befasst, 
eine Fille wertvollster Anregungen bietet. H. J. BUrcGin 


Integral Equations and Their Applications to Certain Problems in 
Mechanics, Mathematical Physics and Technology. Von S.G. MIKuHiin. 
Ubersetzt aus dem Russischen von A. H. ARmstronG (Pergamon Press, London, 
New York und Paris 1957). 338 S., 24 Fig.; 80s. 

Diese Ubersetzung zeigt erneut das hohe Niveau der russischen Lehrbiicher, 
wo vor allem dem methodischen Aufbau der grésste Wert beigemessen wird. Es 
geht dem Verfasser darum, das enorme Material, das sich im Verlaufe der vergan- 
genen dreissig Jahre in Fachzeitschriften angehauft hat und welches sich mit den 
Anwendungen von Integralgleichungen in Physik und Technik befasst, zu sammeln 
und in einer derartigen Form zusammenzufassen, dass es getrost als ein Standard- 
werk auf diesem Gebiet angesehen werden kann. 

Der erste Teil des Buches befasst sich mit den fundamentalen Theorien iiber © 
Integralgleichungen und mit den Methoden der approximativen Lésungsverfahren. 
Spezielle Beachtung wird den singularen Integralgleichungen gewahrt, welche wohl 
schon griindlich erforscht und bei vielfachen Anwendungen ausgepriift sind, jedoch 
noch kaum in einem systematischen Aufbau dargestellt wurden. Wenn immer 
moglich werden die theoretischen Resultate durch numerische Beispiele erlautert. 
Der zweite, gréssere Teil ist den Anwendungen der Integralgleichungen gewid- 
met, wobei hauptsachlich ein- und zweidimensionale Probleme aus den Gebieten der 
Elastizititslehre (biharmonische Probleme) und der Hydrodynamik behandelt 
werden. Ein spezielles Kapitel befasst sich mit den Anwendungen der Theorie der 
singularen Integralgleichungen. TH. GINSBURG 
 Ballistik. Von HERMANN ATHEN, 2. Aufl. (Quelle & Meyer, Heidelberg 1958). 
258 S., 59 Fig., 4 Zahlentafeln; DM 29.-. 

_ Das den Ballistikern wohlbekannte Lehrbuch von H. ATHEN erscheint nun in 
einer zweiten Auflage, die den bedeutenden Entwicklungen wahrend und seit dem 
Kriege weitgehend Rechnung tragt. Die Hauptteile des Buches sind: Aussere 
Ballistik, Treffwahrscheinlichkeit, innere Ballistik, Bomben- und Raketenballistik. 
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Bei der Aussenballistik werden beim Problem des Luftwiderstandes die Grundlagen 
der Strémungslehre und der Gasdynamik gegeniiber der ersten Auflage bedeutend 
ausfiihrlicher behandelt, wahrend die analytischen Lésungsmethoden der (gena- 
herten) ballistischen Hauptgleichungen etwas in den Hintergrund getreten sind, 
obwohl auch jetzt noch Naherungsmethoden und Naherungsansatze bei verschie- 
denen Fragestellungen einen ziemlich breiten Raum einnehmen. : 

Das Kapitel iiber die Stérungstheorie wiirde man vielleicht gerne noch etwas | 
erweitert sehen durch eine ausfiihrlichere Behandlung der praktisch wichtigen 
hdhenveranderlichen Stérungen, insbesondere auch im Falle der Flak-Bahnen, da 
hier die Verhaltnisse haufig wesentlich anders liegen als bei den im Mindungs- 
horizont endenden Flugbahnen. Die Bomben- und Raketenballistik werden dem 
einfiihrenden Charakter des Buches gemdss nur kurz gestreift. Bas Buch bietet 
eine recht brauchbare Einfiihrung in die Fragestellungen der Ballistik, wobei die 
theoretischen Entwicklungen im Vordergrund stehen, wahrend die experimentellen’ 
und messtechnischen Probleme meist nur knapp angedeutet werden. Uber 250 Li- 
teraturangaben erméglichen es aber dem Interessierten, auf ausfiihrlichere Spezial- 
untersuchungen zuriickzugreifen. Das Buch ist ohne Zweifel — wie schon die erste 
Auflage — gut geeignet, manchem die ersten Grundlagen der theoretischen Ballistik 
zu vermitteln. E. RotH-DESMEULES” 


" 
Anfangswertprobleme bei partiellen Differentialgleichungen. Von 
R. SAUER (Springer-Verlag, Berlin 1958). 284 S., 68 Abb.; DM 38.-. : 
Dieses in der «gelben Reihe» des Springer-Verlags (Grundlehren der mathema- 
tischen Wissenschaften) herausgekommene Standardwerk erscheint in der zweiten 
Auflage mit einer Reihe von Erganzungen, welche zum Teil auf den neuesten Ent- 
wicklungen in der Mathematik beruhen. Dies gilt vor allem von dem neu aufge- 
nommienen Kapitel iiber die Behandlung von Anfangswertproblemen mit Hilfe des. 
Distributionskalkiils von LAURENT SCHWARTZ, welcher erst in den vergangenen zehn 
Jahren entwickelt wurde. 
Nach einer kurzen Darlegung der Grundbegriffe der Distributionstheorie, wobei 
unter anderem auch die Pseudofunktionen von H. Riesz eingefiihrt werden, wird fe 
Kalkiil auf Faltungsgleichungen, und zwar besonders auf Anfangswertprobleme der 
Wellengleichung angewandt. In einem weiteren Abschnitt wird die Brauchbarkeit 
des Kalkiils an praktischen Beispielen aus der Theorie der dreidimensionalen lineari- 
sierten Uberschallstr6mung um flache und schlanke Flugkérper erlautert. Hierbei 
erscheint die Hadamardsche Theorie der endlichen Bestandteile von Integralen als 
blosser Spezialfall der Distributionstheorie. 
Bei der Erweiterung der Laplace-Transformation auf Distributionen bleiben die 
Rechenregeln der Laplace-Transformation im wesentlichen giltig, jedoch treten 
gewisse Vereinfachungen ein, und verschiedene einschrankende Voraussetzunge 
fallen weg. So ist es nach der Erweiterung des Oberraumes durch Distributionen nun 
durchaus modglich, die Laplace-Transformation auch auf Dirac-Distributione 
(Sprungfunktionen) anzuwenden. 
Das bereits in erster Auflage bewahrte Buch hat durch diese Erganzungen ein 
wesentliche Bereicherung erfahren. TH. GINSBUR 


Einfuhrung in die Physik. Band II: Elektrizitdtslehve, Wellenlehve, Akusti 
und Optik. Von P, FRAUENFELDER und P, HuBER (Ernst Reinhardt Verlag Miin 
chen und Basel 1958), 504 S., 697 Fig.; Fr. 29.50. 

Der mit Spannung erwartete zweite Band der Einfiihrung in die Physik voi 
FRAUENFELDER und HUBER umfasst die Elektrizitatslehre und die Wellenlehre mit 
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Akustik und Optik. Wiederum iiberrascht die reiche Fiille des Stoffes und die wohl- 
dokumentierte und verstandliche Darstellung mit den vielen, die Anschaulichkeit 
und das Verstandnis férdernden Beispielen. Der Fachmann wird die Moglichkeit 
begriissen, Altbekanntem in neuen Formen zu begegnen, der Lernende aber, fiir 
den das Buch in erster Linie gedacht ist, wird die klare, ausfiihrliche und sich doch 
immer auf das Notwendige beschrankende Darstellung zu schatzen wissen. 

Uber Einzelheiten der Darbietung des Stoffes mag man in guten Treuen anderer 
Auffassung sein, in manchen Abschnitten hatte eine sorgfaltige Ausfeilung vielleicht 
noch zu einer besseren Konzentration des Stoffes fiihren kénnen, und unbefriedigend 
geblieben ist wie im ersten Band die typographische Gestaltung der Formeln. Man 
wird es ferner bedauern, dass in einem: 1958 erschienenen Buche, in welchem dem 
Mechanismus der elektrischen Leitung 70 Seiten reserviert sind, die Halbieiter- 
eigenschaften iiberhaupt nicht erwahnt werden, bedauern wird man auch, dass 
neben dem abstrakten Begriff der Ladungsquelle Angaben iiber Einzelheiten der 
Stromerzeugung vollig fehlen und damit auch der immer noch wertvolle Akkumu- 
lator nirgends Erwahnung findet. 

Durch diese kleinen Einwande soll aber der ausgezeichnete Gesamteindruck des 
ganzen Buches in keiner Weise eine Einschrankung erfahren. Jeder Abschnitt 
er6ffnet bei eingehendem Studium neue, interessante Gesichtspunkte und zeugt fiir 
die manchmal vielleicht etwas eigenwillige, aber immer originelle Darstellungsart 
der beiden Verfasser. Auf besondere Liebe zur Sache und besondere Sachkenntnis 
weisen vor allem manche in den letzten Einzelheiten sehr sorgfaltig redigierte 
Kapitel der Optik hin. 

Die grosse und iiberaus verdankenswerte Leistung, die in der Umarbeitung der 
umfangreichen Stoffgebiete der heutigen Physik in eine fiir Studierende didaktisch 
brauchbare Form liegt, kann der Aussenstehende wohl kaum voll ermessen — ein 
nur dusserliches Zeichen dieser Leistung ist schon die Zahl von beinahe 700 Abbil- 
dungen des vorliegenden Bandes. Man wird daher Verstandnis dafiir aufbringen 
miissen, wenn der mit grossem Interesse erwartete 3. Band, der die Atomphysik 


enthalten wird, nicht schon im laufenden Jahr herausgegeben werden kann. 
W. HARDMEIER 


Ordinary Difference-Differential Equations. Von EpmuND PINNEY 
(University of California Press, Berkeley 1958). 262 S., 45 Fig.; $5.00. ‘ 

Eine Differenzen-Differentialgleichung besteht aus einem Ausdruck, welcher 
sowohl eine Funktion y,, 4, ...4;,(¢) und einige ihrer Ableitungen beziiglich t als auch 
Differenzen beziiglich den h’s enthalt. Ein Beispiel einer solchen Gleichung ist 


Vp lt) = Vn—r(t) - 

Das vorliegende Buch befasst sich mit gewohnlichen (im Gegensatz zu partiellen) 
Differenzen-Differentialgleichungen. Nach einem einfiihrenden und den Problem- 
kreis beschreibenden Kapitel entwickelt der Verfasser die zugehérige Theorie, 
welche anhand von numerischen Beispielen erlautert wird. TH. GINSBURG 
Max-Planck-Festschrift 1958. Herausgegeben von B. KockEL, W. MACKE 
und A. Papapetrou. (VEB Deutscher Verlag der Wissenschaften, Berlin 1959). 
mi2 S:, 28 Abb.; DM 54.-. 

- Die vorliegende Festschrift ehrt Max PLANCK zu Anlass seines hundertsten 
Geburtstags. Sie enthalt 31 Beitrage von 33 Autoren. Naturlicherweise uberwiegen 
die Beitrage aus Europa (25) an Zahl die aussereuropaischen. 
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Die Arbeiten sind nach Inhalt geordnet, wobei neben den Abschnitten iiber_ 
klassische Physik und Relativitaétstheorie vor allem die Quantentheorie mit ihren 
erkenntnis-theoretisch-philosophischen Ausstrahlungen vertreten ist. Die letzteren 
nehmen einen iiberraschend breiten Raum ein (etwa 70 Seiten) und machen auf den 
Rezensenten den ungewohntesten Eindruck. Es ist viel die Rede von der «ortho-— 
doxen Schule der Quantenmechanik», und gegen den Positivismus und den grossen 
Ernst Macn (dessen klassische Mechanik in threr Entwicklung hoffentlich immer 
besinnliche Leser finden wird) wird polemisiert. Man kann sich auch des Eindruckes 
kaum erwehren, dass der Jubilar selbst in einen bekannten weltanschaulichen 
Karren eingespannt werden soll. F 

Sehr beachtenswert ist in diesem Zusammenhang der Artikel von V. Fock 
(Leningrad) Uber die Deutung der Quantenmechanik. In ihm werden mit grosser 
Klarheit von einem wohldefinierten philosophischen Standpunkt aus die grund- 
sdtzlichen Fragen der Quantentheorie durchleuchtet. Das fiihrt zu einer wesentlichen” 
Erweiterung des urspriinglichen Ideenkreises, allerdings immer innerhalb des zu- 
grunde gelegten philosophischen Rahmens. Dass nicht jeder diesen Rahmen aner- 
kennen kann, fallt in Anbetracht der Offenheit der Auseinandersetzung weniger 
ins Gewicht. 

Natiirlich enthalt die Festschrift eine grosse Reihe rein wissenschaftlicher Bei-" 
trage von zum Teil hoher, zum Teil unterschiedlicher Qualitat, auf die der Rezen- 


sent nicht weiter eingehen kann. REs Jost 


gelung. Von A. I. Lurye. Ubersetzt aus dem Russischen von H. KINDLER und 
R. Reissic (Akademie-Verlag, Berlin 1957). 167 S., 24 Abb.; DM 15.-. 

In diesem Buche sind einige der Arbeiten des Verfassers zusammengefasst, die. 
in den Jahren 1945-1950 in russischen Zeitschriften erschienen sind. . 

Im ersten Kapitel wird die Problemsteillung gegeben und eine Methode ent- 
wickelt zur Reduktion der Bewegungsgleichungen eines selbsttatig geregelten 
Systems auf eine Form, die der Verfasser als kanonisch bezeichnet. Im zweiten 
Kapitel werden die Kriterien fiir die Stabilitat «im Grossen» von nichtlinearen 
Regelungssystemen des betrachteten Typs aufgestellt; das dritte Kapitel enthalt 
eine Untersuchung der Eigenschwingungen in diesem System, wahrend sich das 
vierte Kapitel mit dem Verhalten eines Systems an der Grenze seiner linearen 
Stabilitat beschaftigt. 

Das Buch wendet sich in erster Linie an den mathematisch gebildeten Ingenieur, 
der sich mit Berechnung und Bau von selbsttatigen Regelungssystemen befasst, 
aber auch an Physiker und Mathematiker. TE. GINSBURG 


; 
Einige nichtlineare Probleme aus der Theorie der selbsttatigen Re- 
; 


Theoretical Elasticity. Von Cart E. Pearson (Harvard University Press 
Cambridge, Mass. 1959). 218 S., 19 Fig.; $6.-. : 

In der Reihe «Harvard Monographs in Applied Science» ist als dritter Band 
bereits die ausgezeichnete, wenn auch leider Fragment gebliebene Elastizitats- und 
Plastizitatstheorie von H. M. WEsTERGAARD erschienen. Die Herausgabe des vor- 
liegenden Buches als sechster Band der gleichen Reihe stellte daher ein gewisse 
Wagnis dar, das aber in jeder Beziehung gegliickt ist. 

Annahernd die Halfte des Werkes befasst sich mit der Entwicklung der Tensor- 
rechnung unter Beschrankung auf kartesische Koordinatensysteme sowie mit der 
entsprechenden Darstellung der Grundlagen der Kontinuumsmechanik, namlich des 
Spannungs- und des Verzerrungszustandes. Dabei beginnt die Entwicklung wie 
alle amerikanischen Biicher, mit den Regeln der Vektoralgebra; sie eee ae 


: 
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Beispiel in den Abschnitten iiber Differentialgeometrie oder Kreiselung, einiges 
was uber das unmittelbare Ziel des Buches hinausgeht. In der zweiten Halfte 
werden die Grundgleichungen der Elastizitatstheorie formuliert und alle grund- 
satzlich wichtigen Verfahren zur Gewinnung von Lésungen besprochen, wobei ein 
besonderes Kapitel den Variationsprinzipien gewidmet ist. Die drei letzten Kapitel 
behandeln Thermoelastizitat, elastische Wellen und nichtlineare Elastizitat. 

Die zweite Halfte dieses Buches demonstriert eindriicklich die Okonomie der 
modernen Tensorschreibweise. Es ist hier auf 136 Seiten eine erstaunliche Fiille von 
Stoff vereinigt. Zwar werden die bei ebenen Problemen heute im Vordergrund 
stehenden funktionentheoretischen Verfahren nicht dargestellt, und zudem fehlen 
in diesem zweiten Teil (im Gegensatz zum ersten) alle Anwendungen sowie Ubungs- 
aufgaben. Gerade dadurch ist aber auf wenig mehr als 100 Seiten ein ausserst niitz- 
liches und iibersichtliches Kompendium der raumlichen Elastizitatstheorie ent- 
standen. H. ZIEGLER 


Lehrbuch der Hochfrequenztechnik. Von Fritz Vivsia, 5. Aufl., Band II 
(Akademische Verlagsgesellschaft mbH, Frankfurt am Main 1959). 735 S., 950 Abb.; 
DM 48.-. 

Das Lehrbuch der Hochfrequenztechnik von ViLBIG, das 1937 seine erste Auflage 
erlebte, kommt nun in fiinfter, neubearbeiteter Auflage heraus; diese umfasst drei 
Bande, von denen der mittlere 1958 erschienen ist. Das Lehrbuch hat sich zu einem 
eigentlichen Standardwerk entwickelt, das dank seiner Vollstandigkeit und Uber- 
sichtlichkeit den Dienst eines Handbuches versieht und der Fachwelt unentbehrlich 
ist. Von der vorbildlichen Qualitat zeugt der hier zu besprechende zweite Band der 
neuesten Auflage. Er umfasst in reicher Fiille: Verstarker (inkl. réhrenlose Verstar- 
ker und Laufzeitverstarker), Schwingungserzeugung, Riickkopplung (inkl. Gegen- 
kopplung und Neutralisation), Abschirmung, Gleichrichtung, Modulation und 
Demodulation, Tastung und Impulsmodulationsverfahren, Frequenzwandlung, 
Schwundregelung, Dynamikanderung sowie einen Nachtrag (molekulare Mikro- 
wellenverstarker und Oszillatoren). Die Ausstattung ist, wie bei den friiheren Auf- 
lagen, sehr gut; ausserordentlich zu schatzen sind die zahlreichen Abbildungen. 

F. Tank 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von — 
S. FLtecr. Band 12: Thermodynamik dey Gase (Springer-Verlag, Berlin 1958). 
686 S., 253 Fig.; DM 154.-. 

Im Rahmen des durch S. FLiGGE neu herausgegebenen Handbuches der Phystk 
ist der Thermodynamik der Gase ein besonderer Band gewidmet. Darin werden die 
verschiedenen Aspekte dieses Wissensgebietes in umfassender Weise behandelt. Die 
verschiedenen Hauptkapitel und ihre Verfasser sind die folgenden (die Beitrage 
sind teils in deutscher, teils in englischer Sprache abgefasst) : 

The Properties of Real Gases (J. S. RowLtnson) ; 

Theory of Real Gases (J. E. Maver); 

Principles of the Kinetic Theory of Gases (H. GRAD); 

Transporterscheinungen in Gasen von mittlerem Druck (L. WALDMANN) ; 
Allgemeine Vakuumphysik (R. JAECKEL); 

Production and Measurement of Ultrahigh Vacuum (D. ALPERT). 

Z Vorweg sei bemerkt, dass es wohl kaum méglich sein wird, einen Referenten zu 


oa 


finden, der zur Beurteilung aller Teile dieses Buches in gleicher Weise zustandig 
ware. Hingegen zeigt schon jede Stichprobe, dass dieses Nachschlagewerk ein 
ausserordentlich hohes wissenschaftliches Niveau halt. 

i4 
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Der Beitrag von RowLINson tiber die Eigenschaften der realen Gase halt sich 
im Rahmen der klassischen phanomenologischen Thermodynamik und beleuchte 
sowohl die theoretischen Formulierungen (Zustandsgleichungen usw.) als auch die 
experimentellen Methoden und empirischen Ergebnisse. Eine besonders sorgfaltig 
Behandlung erfahrt das Gebiet des kritischen Punktes, auch fiir den Fall der Gas-_ 
mischungen. 

Der von Maver verfasste Teil bringt die theoretischen Grundlagen von hoher 
wissenschaftlicher Warte aus zur Darstellung, wobei selbstredend vom statistische 
(auch quantenstatischen) Standpunkt ausgegangen wird. Im vierten Hauptabschnitt 
dieses Teiles wird auf das Quantengas eingegangen. ’ 

Die kinetische Theorie selbst, die bekanntlich vor allem mit dem Namen BoLtz - 
MANNS verkniipft ist, erfahrt eine ausfiihrliche und tief in die begrifflichen Grund- | 
lagen eindringende Behandlung durch GRAD. 

Der Artikel von WALDMANN iiber die Transporterscheinungen in Gasen (Zahig- 
keit, Warmeleitfahigkeit, Diffusion, Thermodiffusion) nimmt entsprechend der 
hohen praktischen Bedeutung dieses Gebietes etwa einen Drittel des ganzen Buches 
ein. Die phanomenologische Betrachtungsweise wird zum Ausgangspunkt genom= 
men, doch wird spater die tiefere Begriindung aus der kinetischen Theorie gewon- 
nen. Experimentelle Methoden und Ergebnisse erfahren eine ausfiihrliche Wiir- 
digung. 

Die letzten zwei Hauptteile des Buches schliesslich befassen sich mit der Vakuum- 
physik. Sie unterscheiden sich von den vorhergehenden namentlich dadurch sehr 
deutlich, dass hier der Schwerpunkt auf der Behandlung der technischen Vorrich 
tungen (Pumpen, Messgerate) liegt. 

Jeder der Verfasser hat das von ihm bearbeitete Gebiet als ein geschlossenes 
Ganzes behandelt. Dabei werden natiirlich gewisse Uberschneidungen unvermeid- 
bar, was aber nicht als ein Nachteil des Buches gewertet werden darf. Auch der 
durchaus individuelle Charakter der einzelnen Beitrage ist als Positivum zu betrach- - 
ten. Alles in allem entspricht das Buch den hohen Erwartungen, mit denen man es 
von vornherein zur Hand nimmt. Es wird ein hervorragendes Nachschlagewe 
sein fiir einen weiten Kreis von Forschern. Zu diesen gehéren vor allem Physik 
und Physikochemiker, im Zuge der heutigen Entwicklung aber bereits auch Inge= 
nieure, was in Zukunft noch in h6herem Masse der Fall sein wird. W. TRAUPEL 


Modern Computing Methods. Notes on Applied Science No.16. Heraus- 
gegeben vom National Physical Laboratory (Her Majesty’s Stationery Office, 
London 1957). 129 S.,18 Fig.; 10s. 6d, 7 

Die unbekannten Verfasser beschreiben in konzentrierter Form numerische 
Methoden fiir Probleme aus verschiedenen Gebieten der angewandten Mathe- 
matik. Beispielsweise wird fiir lineare Gleichungssysteme der GauSsche Alg 
rithmus vorgeschlagen, und es wird (in Kapitel 2) auch gezeigt, wie man d 
Elimination mit einem Rechenautomaten zweckmassig durchfiihrt. Entsprechen¢ 
werden Methoden zur Lésung algebraischer Gleichungen und Eigenwertproble me 
angegeben (Kapitel 3 und 4). Weitere Abschnitte sind der Differenzenrechnu 
und der numerischen Integration von gewéhnlichen sowie partiellen Differential- 
gleichungen gewidmet. Am Schluss findet der Leser noch interessante Ausfiihru ne | 
gen uber die Berechnung und Tabellierung von Funktionen mit Rechenautomaten.. 
Wenn auch neuere Methoden wie konjugierte Gradienten, Bi-Iteration usw. 
fehlen, gibt das kleine Werk doch wertvolle Hinweise fiir dea Mathematiker 
der mit Rechenautomaten zu tun hat. H RurisHausHl 
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